The Annals of Probability 

2007, Vol. 35, No. 6, 2160-2212 

DOI: 10.1214/009117907000000141 

© Institute of Mathematical Statistics, 2007 

SECOND ORDER ASYMPTOTICS FOR MATRIX MODELS 

By Alice Guionnet and Edouard Maurel-Segala 

Ecole Normale Superieure de Lyon 

We study several-matrix models and show that when the poten- 
tial is convex and a small perturbation of the Gaussian potential, the 
first order correction to the free energy can be expressed as a gener- 
ating function for the enumeration of maps of genus one. In order to 
do that, we prove a central limit theorem for traces of words of the 
weakly interacting random matrices defined by these matrix models 
and show that the variance is a generating function for the number 
of planar maps with two vertices with prescribed colored edges. 

1. Introduction. In this paper we study the asymptotics of Hermitian 
random matrices whose distribution is given by a small convex perturbation 
of the Gaussian Unitary Ensemble (denoted GUE). We shall consider m- 
tuples of random matrices, with an integer number m € N fixed throughout 
this paper. Then, the law /x of m independent matrices following the GUE 
is given, for N x N Hermitian matrices A = (A±, . . . , A m ), by 

m N 

dfj, N (A) = e (-W*(E£i^) JJ JJ diAAjj J] dSte^AA^dShniAA^, 

i=lj=l l<j<k<N 

with Tr the nonnormalized trace Tr(A) = J2iLiAu- In other words, the 
A = (A\, . . . ,A m ) are independent Hermitian matrices whose entries are, 
above the diagonal, independent complex centered Gaussian variables with 
variance A -1 . Let V(X) be a polynomial in m noncommutative indeter- 
minates X = (X±,..., X m ) such that Tr(V(A)) is real for all m-tuple of 
Hermitian matrices A = (A±, . . . , A m ). Then, we shall study the following 
probability measure jiy on the set WAr(C) m of m-tuple of N x N Hermitian 
matrices 

d^(A) = -L e-^(A) )V v (A)) 
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where Zy is the normalizing constant so that fj,y is a probability measure. 

Besides, we require that the trace of W(A) := V(A) + \ YaLi Af is a 
strictly convex function of the entries of A = (A±, . . . ,A m ) E T-L^{C) m for 
any N E N. In that case Zy is automatically finite. More precisely, for c > 0, 
we say that V is c-convex if for any N E N, A E W7v(C) m -> Tr(W(A)) is 
real-valued and with Hessian bounded below by ci\ An example of c-convex 
potential is 

V(X 1 , . . . , X m ) = Pj ( E 4^ J + E ^'■■.V;A- ; , 

j \ i ) j,k 

with convex polynomials Pj on M, real numbers ct\ , /^fc and J2j \Pj,k | < 1 — c 
for all k E {1, . . . , m} (see Section 2 for more details). 

The central result of this paper can roughly be stated as follows. 

Theorem 1.1. Let V = V t (X 1 ,...,X m ) = £"=i t jgj (X 1} ...,X m ) be a 
polynomial potential with n G N, t = (t\, . . . , t n ) S C n and monomials (qj)i<j< 
fixed. For all c > 0, i/iere exists rj > so £/ia£ i/ |t| := maxi<j< n \ tj\ <rj and 
Vt is c-convex, there exists F l (Vt) = F l (t\, . . . , t n ) for i = 0, 1 so £/iaf 

logZ^ = iV 2 F°(y t ) + F 1 ^) + o(l). 

The first order expansion F°(Vt) was already obtained in [14] and we 
extend our study here to the second order. The higher order expansions 
can also be tackled by a refinement of our strategy; this is the subject of 
a separate article by Maurel-Segala [19]. Moreover, we believe our tools 
sufficiently robust to tackle other models such as the Gaussian orthogonal 
ensemble, or the Haar measure on the unitary group, for instance. Again, 
this is the subject of further studies. 

We next turn to the combinatorial interpretation of F°(Vt), i ?1 (Vt) has 
generating functions of maps. 

Matrix models have been used intensively in physics in connection with 
the problem of enumerating maps; see the reviews [9, 12]. Let us recall that 
a map of genus g is a graph which is embedded into a surface of genus g in 
such a way that the edges do not intersect and dissecting the surface along 
the edges decomposes it into faces which are homeomorphic to a disk. We 
will call a star the couple of a vertex and the half-edges which are glued to 
this star. A star will have a distinguished half-edge and an orientation and 
will eventually have colored half-edges when m > 2. When m = 1, it is well 
known that if Vt = [i.e., t = (0, ... ,0)] moments of the random matrices 
from the GUE are related with the enumeration of maps; for instance, the 
number M 9 k of maps with genus g with one star with 2k half-edges were 
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computed by Harer and Zagier [16] using the formula 




It was shown in [10] (see also [1, 2]) that when m = 1, this enumerative prop- 
erty extends to the free energy of matrix models at all orders, as conjectured 
and widely used in physics (see, e.g., [7]). More precisely, if Vt = Y^i=iti% ni 
with D = maxHi = n p even and t p / \U\ large enough, for all k G N, there 
exists r\ > so that, for |t| < i], 

log^ = iV 2 E^(^t) + (iV 2 - 2fc ), 

3=0 

with 

*"W = E W^mi^ 

fci,...,fc„eN"\{0,...,0} i l ' 

where M 9 ki kn is the number of maps of genus g with ki vertices of degree 
n i, 1 <i <n. 

Several-matrices integrals are related with the enumeration of colored 
(or decorated) maps. To make this statement clear, let us associate to a 
monomial q(X) = X{ 1 ■ ■ ■ Xi p a colored star as follows. We choose m different 
colors {1, . . . , m}. The star associated to q (called a star of type q) is a vertex 
equipped with colored half-edges such that the first half-edge has color i\, 
the second has color ii till the last half-edge which has color i p . Because 
the star has a distinguished half-edge (the one associated with X^) and 
an orientation, this defines a bijection between noncommutative monomials 
and colored stars. Then, it can be seen [23] that, for any monomial q, 

hm / ±-Tr(q(A))dn N (A)=M (q), 

with A4o(q) the number of planar maps with one-colored star of type q such 
that only half-edges of the same color can be glued pair-wise together (then 
forming a one-colored edge). In [14], we proved that if Vt is c-convex and 
t = (<!,.. . ,t n ) is small enough, the limit F°(Vt) of the free energy given in 
Theorem 1.1 is analytic in the variables ti in a neighborhood of the origin 
and its expansion is a generating function for planar maps with prescribed 
colored stars; 

(i) F°(v t )= E n^r^,..,^ 

fei,...,fe„6N"\{0,...,0}i=l 
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with Adk lt ...,k n the number of planar maps with ki colored stars of type q^ 
the gluing being allowed only between half-edges of the same color. Note, 
however, that we cannot retrieve all the numbers Mk\,...,k n from the i ?0 (Vt)'s 
because the condition that Tr(T4) is real requires that the parameters t 
satisfy some relations. Namely, if * denotes the involution (zX^ ■ ■ ■ X{ k )* = 
zXi k we must have Tr(Vt) = ^Tr(Vt + Vj*) and, therefore, if 14 = 

hO.ii to each must corresponds a tj such that Tr(g^) = Tr(g*) and tj = ti- 
Thus, the F°(Vt)'s are generating functions for the number of planar maps 
with ki colored stars of type qi or q* . The convexity assumption also should 
induce some extra relations between the parameters, but it can be removed 
as shown in Theorem 1.4. 

In this paper we shall prove that such a representation also holds for the 
correction F 1 (Vt) to the free energy given in Theorem 1.1. 

Property 1.2. F 1 (Vt) is analytic in the parameters ti in some neigh- 
borhood of the origin. Its expansion is a generating function of maps: 

e n^?X,..,* n , 

fci,...,fc n eN™\{0,...,0}i=l l ' 

with M\ k n the number of maps with genus one with ki colored stars of 
type q$. In particular, the above sum converges absolutely for maxi<j< n |tj| 
small enough. 

Let us remark that such a representation is commonly assumed to hold 
in physics since the formal result is always true for finite N. For a few 
models (viz., models similar to the Ising model on random graphs), the 
analysis has been pushed forward to actually give a rather explicit formula 
for the generating function i ?1 (\4) m terms of the limiting spectral measure 
of one matrix under the Gibbs measure fiy (see, e.g., Eynard, Kokotov and 
Korotkin [11]). Our strategy is here to study the most general potentials, 
providing a general formula for i ?1 (Vt) in terms of the limiting empirical 
measure of all the matrices (see Section 6). 

Our arguments to prove Theorem 1.1 are rather different from [10] or [1] 
where orthogonal polynomials were used. In [10], the idea was to develop a 
Riemann-Hilbert approach based on precise asymptotics of orthogonal poly- 
nomials. In the case of several-matrices models, the technology of orthogonal 
polynomials is far to be as much developed (except for the Ising model; see 
[6]). We shall therefore use different tools; the first, which is well spread in 
physics, is the use of the Schwinger-Dyson equation, the second, for which 
we need a convex potential, is the a priori concentration inequalities. To 
sketch our strategy, let us denote p, the empirical measure 

V N : P — Itt(P(A)) = ±Tr(P(A u A m )), 
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where P runs over the set C{X\, . . . ,X m ) of noncommutative polynomials 
in m indeterminates. Note that when m = 1, fi N is the spectral measure of 
Ai, and therefore a probability measure on R. When m > 2, jl N is a tracial 
state, which generalizes the notion of measures to a noncommutative setting 
(see, e.g., [24]). Observe that, for 1 < i < m, 

d tt logZ» = -N 2 ^ N ( qi )) 

so that the second order asymptotics of the free energy will follow from 
that of ~p N = fiy [fi N ] evaluated at the monomials qi, 1 < i < n. Then, a 
simple integration by parts shows that, for any ^GN, the following finite 
N Schwinger-Dyson equation holds 

vv t (A* ® A*(W) = Mh(A jv ((^ + AH)P)) 

for any polynomial P and i G {1, .. . , m}. Here, di,D{ are noncommutative 
derivatives (see Section 2 for a definition) . Based on this equation and con- 
centration inequalities, it was shown in [14] that for sufficiently small pa- 
rameters t = (ti, . . . , t n ), fi N converges almost surely and in expectation (for 
the weak topology generated by the set C(Xi, . . . ,X m ) of noncommutative 
polynomials). Its limit fit is a solution of the Schwinger-Dyson equation 

(2) Mt ® fH(diP) = fH{(Xi + DiV t )P) 

\/PeC(x 1 ,...,x m ),i<i<m. 

It is the unique solution which satisfies a bound of the form |/xt(X^)| < C d 
for all d € N and alH G {1, . . . , m}, when t = (t\, . . . , t n ) is small enough and 
C finite, independent of t. 

In this paper we investigate the correction to this convergence by proving 
a central limit theorem for fi N — fif More precisely if we define 5^(P) := 
N(jl N (P) — fit(P)), then we show the following: 

Theorem 1.3. For all c > 0, there exists rj > such that for all t in 
B^ c = B(0,rj) n {t|Vt is c-convex}, for all P in C(X\, . . . ,X m ), under fj,y , 
§t(P) converges in law toward a complex centered Gaussian law jp. More- 
over, {"fp\P G C(Xi, . . . , X m )}, equipped with the natural addition Jp + Jq = 
jp + q, is a Gaussian space and the covariance function is a generating func- 
tion for planar maps with two prescribed stars. 

Such a central limit theorem was proved for more general potentials when 
m = 1 by Johansson in [18]. When m = 1 but the entries are not Gaussian, 
we refer the reader to [3]. In the case m>2 but V = 0, the central limit 
theorem was obtained in [8, 13, 20]. Our proof is rather close to that of [18] 
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and in the physics spirit; by doing an infinitesimal change of variables, it 
can be seen that the random variable 

m 

$t(&P) ■= E $t((I ® Mt + Mt ® I){diDiP) - (X t + AH)AP) 

i=l 

converges in law toward a centered Gaussian variable. The main issue is then 
to show that the £t-f"s are dense in the set of polynomials. When m = 1, 
Johansson could use finite Hilbert transformation to invert the operator £f 
In our case, we deal with a differential operator acting on noncommutative 
test functions and we prove by hand that it is invertible for sufficiently small 
tj's in Section 4. Clearly, our analysis is perturbative at this point and does 
not try to find the optimal domain of validity of the central limit theorem. 

To use the central limit theorem to obtain the second order asymptotics of 
fj,y (5t(P)) observe that, by the finite dimensional Schwinger-Dyson equa- 
tion, we get 

Nf*v t 0?foP))=$ t (S?®8?(a i D i P)) 

and the right-hand side converges toward the variance of the central limit 
theorem. So again, to obtain the limit of Ny$ x (8?(P)), we need to invert the 
operator £t ( see Section 6). The resulting formula for the free energy and the 
variance are given in terms of differential operators acting on noncommuta- 
tive polynomial functions. Note that a similar formula for the variance of the 
central limit theorem governing the fluctuations of words of band matrices 
was found in [13] . Their interpretation in terms of enumeration of maps can 
be retrieved from the interpretation of noncommutative derivatives in terms 
of natural operations on maps (see [14]). 

Finally, in the spirit of [14], we study matrix models with a nonnecessarily 
convex potential V. Since in that case Zy has no reason to be finite, we need 
to add a cut-off. For a positive constant L, we define 

^ L (dA) = -^l w(A)<i e-™(vW) dfl N (Ah 

V,L 

with A max (A) the maximum of the spectral radius of the AiS and Zy L 
a normalizing constant. The remarkable point that we shall prove is that 
asymptotically the behavior of this measure is independent of L and gives 
the same type of expansion as in the convex case. 

Theorem 1.4. Let V = V t {X 1} ..., X m ) = £? =1 t jQj (X u . . . , X m ) be a 
polynomial potential with n € N, t = (t±, . . . , t n ) G C n and monomials (qj)i<j< 
fixed. Assume that Tr (14(A))) is real for all A G H N (C) m , all N G N. There 
exists Lq > such that, for all L > Lq, there exists rj > so that if |t| := 
maxi<j< n \tj \ <n then 

log Z$ t>L = N 2 F°(V t ) + F\V t ) + o(l) 
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with F° (Vt) , F 1 (Vt) as in (1) and Property 1.2. 

Note that in the large N limit, the dependence in L disappears. 

In the next section we will describe our hypothesis of convexity and show 
some useful consequences. In Section 3 we give an estimate on the rate 
of convergence of fj,y [fi N ] to /it . Then, in Section 4 we prove a central 
limit theorem, first only for some specific polynomials and then for arbitrary 
polynomials. In Sections 5 and 6, we give an interpretation of the variance 
and of the free energy in terms of enumeration of maps. Finally, in Section 
7 we give some hints to generalize our proofs to the setting of Theorem 1.4. 

2. Convex hypothesis and standard consequences. 

2.1. Framework and standard notation. 

2.1.1. Noncommutative polynomials. We denote C{X\, . . . ,X m ) the set 
of complex polynomials on the noncommutative unknown X±, . . . ,X m . Let 
* denote the linear involution such that for all complex z and all monomials 

(zX^ ■ ■ ■ Xi p )* = zXi v ■ ■ ■ Xi x . 

We will say that a polynomial P is self- adjoint if P = P* and denote <C{X\ , . . . , 
X m ) sa the set of self-adjoint elements of C(Xl, . . . ,X m ). 

For an integer number N, we denote TCn(C) the set of N x N Hermitian 
matrices. We shall sometimes identify Wtv(C) with the set W N of the cor- 
responding real entries [by the bijection which associates to A £ 'Hn(c) the 
iV 2 -tuple ((Ke(^)i<i<j<Ar, (^sm(A i j) 1 < i< j< N )))}. 

Moreover, we shall denote in general by A a random matrix, by X a 
generic noncommutative indeterminate (e.g., to write polynomials). Bold 
symbols will in general denotes vectors; A (resp. X) will in general denote a 
m-tuple of matrices (resp. noncommutative indeterminates), whereas t will 
denote a vector of complex scalars. 

The potential V will be later on assumed to be self-adjoint which guar- 
antees that, for all integer JV, all A = (A 1 ,...,A m ) £ H N (C) m , Tr(V(A)) 
is real. Note that, conversely, if Tr(V(A)) is real, Tr(V(A)) = Tr((V + 
V*)(A)/2) and so we can replace V by (V + V*)/2 without loss of gen- 
erality. 

We shall assume also that V satisfies some convexity property in this 
paper. Namely, we will say that V is convex if, for any JV6N, 

<$ : n N (C) m ~ (8 ff2 ) m — ► R, 

((^4fc)ii)i<i<i<iv — ► Tr(V(Ai,. . .,A m )) 

l<k<m 
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is a convex function of its entries. 

Note that as we add a Gaussian potential \ YaLx Xf to V we can relax the 
hypothesis a little. We will say that V is c-convex if c > and V + J2T X? 
is convex. Then the Hessian of 4>w W1 ^h W = V + | J2T Xf is symmetric 
positive with eigenvalues bigger than c. 

An example is 

n / m \ 

v t (x 1 ,...,x m ) = ^pAYl a i x * + E Pk,iX k Xi , 

i=l \k=l ) k,l 

with convex real polynomials Pi in one unknown and real a\,/3k,i such that, 
for all I, J2k \Pk,l\ < (1 — c )- This is due to Klein's lemma (see [15]) which 
states that the trace of a real convex function of a self-adjoint matrix is a 
convex function of the entries of the matrix. 

In the rest of the paper we shall assume that V is c-convex for some 
c> fixed. We will denote B(0, rf) = {t G C n : maxi<j< n |tj| < rf\ and B v>c = 
B(0,rj) H {t : Vt is c-convex}. 

2.1.2. Noncommutative derivatives. We define for 1 < i < m the noncommu- 
tative derivatives di : C(X\,. . . , X m ) — > C(Xi, X m )® 2 by the Leibniz rule 

diPQ = diP x (1 ® Q) + (P ® 1) x diQ 

and c\X,- = li—jl <S> 1. So for a monomial P, the following holds: 

$i> = ^ P g) 5, 

where the sum runs over all possible monomials R, S so that P decomposes 
into RXiS. We can iterate the noncommutative derivatives; for instance, 
df :C(X u ...,X m ) ->C(Xi,...,X m ) ®C(X u ...,X m )®C(X 1 ,...,X m ) is 
given on monomial functions by 

dfP = 2 R®S®Q. 

P=RX i SX l Q 

We denote ft : C{X U . . . , X m )® 2 x C{X U . . . ,X m ) ^ C{X U . . . ,X m ) the map 
P®Q$R = PRQ and generalize this notation to P®Q®R$(S,T) = PSQTR. 
So diPftR corresponds to the derivative of P with respect to Xi in the direc- 
tion R, and, similarly, 2~ 1 [D' 2 P$(R, S) + D?P$(S, R)) the second derivative 
of P with respect to Xi in the directions R, S. 

We also define the so-called cyclic derivative Dj. If m is the map m{A® 
B) = BA, let us define Di = modi. For a monomial P, DiP can be expressed 
as 

DiP = SR - 

P=RXiS 

We shall denote in short D the cyclic gradient (D\, . . . ,D m ). 
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2.1.3. Noncommutative laws. For (Ax, . . . , A m ) € T^at(C)" 1 , we define the 
linear form Aai,..,a to on c ( x i, ■ ■ -, x m) by 

^ u ..., A jP) = ^Tr(P(A 1 ,...,A m )), 

with Tr the standard trace Tr(A) = J2f=i^a- A m w ^ n sometimes be 

called the empirical distribution of the matrices (At, . . . ,A m ). When there 
is no ambiguity and the matrices Ax,..., A m follow the law fXy , we shall 
drop the subscript Ax, ... , A m ; £l n = A^i ... A m - ^ n [^] ^ was shown that if 
T4 = J2itiQi is c-convex, for |t| := maxi<j< n \t{\ small enough, fi N converges 
weakly in expectation and almost surely under fj,y toward a limit fit [i-e., for 
all P in C(Xi, . . . ,X m ), fi N (P) converges in expectation and almost surely 
to Ht(P)]. We denote 

T^(P)=^ N (P)]- 
We shall later estimate differences between p, and its limit. So, we set 

T t N = fFdi# = NZp? -fit), 
If = N(fi» 

In order to simplify the notation, we will make t implicit and drop the 
subscript t in the rest of this paper so that we will denote ~p N , fi,5 N ,8 and 
S_ in place of fl^, fi t , 5^ , 5t N and 6_ t , as well as V in place of Vt- 

2.2. Brascamp-Lieb inequality and a priori controls. We use here a gen- 
eralization of the Brascamp-Lieb inequality shown by Harge [17] which im- 
plies that if V is c-convex, for all convex function g on (M) mAr2 ~ 7i]\f(C) m , 

(3) J g(A-M)d^(A)< J 5 (A)d/if(A), 

where M = J A dfiy(A) is the m-tuple of deterministic matrices with entries 
(Mj)ki = I(Aj)ki dfiy(A) for k,£e{l,...,N}Je{l,.. .,m}, and fi" is the 
Gaussian law on W/v(C) m with covariance (Nc)~ l , that is, / / (A) dfi^f (A) = 
J f(c- 1 / 2 A)dfi N {A) for all measurable function / on W nN2 . 

Recall that P>r\,c is the subset of the complex numbers t E C n which are 
bounded by n and so that V is c-convex. Based on the Brascamp-Lieb 
inequality, the following was shown in [14] (Theorem 3.4): 

Lemma 2.1 (Compact support). Ifc,rj> 0, then there exists Co = Cb(c, 77) 
finite such that, for all i€ {1, . . . , m}, all n £ N, all t £ B„ c , 

»(Xf n )<lim S uplI N (X? n )<Ci n . 

N 
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Note that this lemma shows that, for i G {1, . . . , m}, the spectral measure 
of Xi is asymptotically contained in the compact set [—Co, Co]. 

Proof of Lemma 2.1. Let us recall the proof of this result for com- 
pleteness. Let ke{l,...,m}. As g : A G H N {C) m ■ -» N~ l Ti{Af) = fi N (Xf) 
is convex by Klein's lemma, we can use the Brascamp-Lieb inequalities (3) 
to see that 

(4) H N ((X k - M k f d ) < ^(fi N (Xt d )), 

where M k := Hy{A k ) is the deterministic matrix with entries (M k )ij = fJLy x 
{{Ak)ij). Thus, since fj,^ {fi N [Xf d )) converges by Wigner theorem [25] toward 
c~ 2d C2d < (c _1 4) 2rf with C2d the Catalan number, we only need to control 
M k . First observe that, for all k, the law of A k is invariant under the unitary 
group so that, for all unitary matrices U, 

(5) M k = i4[UA k U*\ = Ui4[A k ]U* => M k = ^(jl N (X k ))I = Ji N (X k )I. 
Let us bound ~fi N (X k ). Jensen's inequality implies 

Z V > e -N^ N ((l/N)Tr(V)) = e -N^ N °V N {V)_ 

According to [23] , fi o fi N converges in moments to the law of m free semi- 
circular operators, which are uniformly bounded. Thus, there exists a finite 
constant L such that > e~ N L . We now use the convexity of V to find 
that, for all N, all A = (A u . . . , A m ) £ H N (C) m , 

(-. m \ / m m \ 

V(A) + Yl A i J > Tr [V(0) + £ DiV(0)Ai + (1 - c) £ Aj . 

By Chebyshev's exponential inequality, and then using the above bound, we 
therefore obtain that, for any A > 0, 

^(fl N (X k )>y)<e-^ 2 y^(e^ N ^) 

rzN 

= -XN 2 y Z c N, XN^ N (X k )-NTr(V(A)+((l-c)/2)yr i A?) ) 
Zj x Ay 



< e N 2 (L-V(0)-\y+m./2\ogc) 

X ( e - JVr& (E™i(( 1 - c )+' D iV(0))^-A^ fc )) 
_ e N 2 (L-V(0)-\y+m/2\ogc) 

x e (N 2 /(2c)) J2^ k (l-c+D l V(0)) 2 +(N 2 /(2c))(l-c+D k V(0)-X) 2 

where we denoted V(0) := V(0, . . . , 0) and DjV(O) := AV(0, . . . , 0). Remark 
that these constants are uniformly bounded for t in B(0,R), R > 0. Thus, 
we deduce that 

/^(/^(A fc ) > y) < e ^[(a-Xy )+ (l/(2c) )( X-b) 2 ] 
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with two constants a, b which are uniformly bounded in terms of c, rj for 
t £ i?r;,c- Optimizing with respect to A shows that there exists A < +00 so 
that, for t in -B^c, 



^(t N (X k )>y)<e N2 («-^y 2 -^ 



,N 2 (A-(c/4)y 2 ) 



Replacing X k by —X k , we bound similarly ^,y(/l N (X k ) < —y) and, hence, 
we have proved 

^(\^ N (X k )\>y)<2e N2 ( A -^y 2 \ 

As a consequence, 

poo 

$(\fi N (X k )\)= ^{\fL N (X k )\>y)dy 

(6) J ° 

<2v^4 + 2 r e-^^-^dy^AV^A, 

where the last inequality holds for iV sufficiently large. Recall that A is a con- 
tinuous function of the U's and, therefore, our bound on sup^r {\fi N (X k )\), 
which controls the spectral radius of M k in any dimension N, is locally 
bounded in t. This completes the proof with (4). □ 

Let us derive some other useful properties due to the convexity hypoth- 
esis. Let A^ ax (^4j) be the maximum of the absolute value of the eigenval- 
ues of A{. We first obtain an estimate on A^ ax (A), the maximum of the 

Ki<m 

under the law fiy. 

Lemma 2.2 (Exponential tail of the largest eigenvalue). Ifc,rj>0, then 
there exists a > and Mq < 00 such that, for all t G B V)C , all M > Mq and 
all integer N , 

^(Al x (A)>M)<e^ AW . 

Proof. Since the largest eigenvalue 

A^ ax (A) = max sup (u,Ai A* u) 1 ^ 2 



KKm 



u =1 



is a convex function of the entries of the Ai's, we can apply the Brascamp- 



— 1 

10 J' 



Lieb inequality (3) to obtain that, for all s G [0, 

3 iVA£UA-M) dfl N {A) < J e sNX^(A) dfl N {A) < ^ 
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where the last inequality comes from the bound on the largest eigenvalue of 
the GUE shown, for instance, in [5]. Now, 

Al x (A) < A^ ax (A - M) + Al x (M) < A^ ax (A - M) + aVIc^, 

where we used the bound (6). Consequently, we deduce that 

Je sNX ^ A U^(A)<C N 

for a positive finite constant C. We conclude by a simple application of 
Chebyshev's inequality. □ 

2.3. Concentration inequalities. We next turn to concentration inequal- 
ities for the trace of polynomials on the subset of H.N(C) m ~ M. N m : 

A M = | A G n N (C) m : A^ ax (A) = max(Al x (^)) < M 

for some fixed M > 0. Recall that 5 = N(jl N — p, N ). We shall prove concen- 

tration inequalities for 5 (P) on A M for polynomial functions P. However, 
concentration inequalities should not restrict to polynomial functions but 
hold more generally for Lipschitz functions (see, e.g., [15]). We thus define 
the following Lipschitz semi-norm: 

/ m \ 1/2 

(7) ||P||f = sup sup [J2\\D k PD k P*\U) ■ 

_4C*-algebra xi,...,x m eA \k=l ) 
Vi,Xi=x* ,\\xi\\^<M 



Be aware that this is not a norm, since, for example, = or [IX1.X2 — 

A^AiH^ = 0. However, on these particular polynomials, <5 vanishes. This 
fact can be generalized as follows; if we set 



m M,P := .,n\n\ Vv(^A% i I N ( P ))> 



we shall see (see the proof below) that on A M 



N 
M 

<T (P) - m^p\ < 2M^/mN\\P\\% . 

Therefore, if we denote C(X±, . . . ,X m )^ the completion and separation of 
C{X\, . . . ,X m ) for || • H^, we can extend <5 — m MP to C(X\, . . . ,X m )j? on 
A similar result will be proved for \i in Lemma 4.9 (note, however, that 
the arguments of this lemma do not apply here because j2 N is not the law 
of uniformly bounded matrices). 
We shall prove the following: 
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Lemma 2.3 (Concentration inequality). Let t he such that V is c-convex. 
There exists a, Mo > such that, for all N in N, all M > Mq, all P G 
C(Xi, . . . ,X m )^ , there exists a positive constant Ep M such that, for any 
£>0, 

(8) ^({|1>) - m N PM \ > e + e N P>M } n A&) < ^Mm??). 
Moreover, there exists a universal constant C such that 

e^ M < 2CNM\\P\\ A c I e- aNM . 
If P is a monomial of degree < d < aN , we have 

||P||f < dM d -\ e^ M < NCdM d e~ aMN , 

\m^ M \<N(3M d + d 2 )e~~ aMN . 

Proof. Since V is c-convex, for all integer number N, the Hessian of 
<$ : A ~ ((^)y)l<i<j<iV G M mJV2 — Tr(y(A l5 . . . , A m )) G R 

l<fc<m 

is bounded below by cl. Therefore, since fiy has density e _Ar< ^v( A ) with 
respect to the Lebesgue measure, fiy satisfies a Log-Sobolev inequality with 
constant (iVc) -1 (see, e.g., Corollaire 5.5.2, page 87 in [4]). In other words, 
for any continuously differentiable function / from M. mN into R, 

with V/ the gradient of / and || • || the Euclidean norm. Here and in the 
sequel we identify the measure fiy as a measure on R^ 2 " 1 . This implies, by 
the well-known Herbst argument (see, e.g., [4], Theoreme 7.4.1, page 123), 
that fiy satisfies concentration inequalities. Let us briefly summarize this 
argument for completeness. If / is a continuously differentiable function, 
differentiating X(X) := j log [ly [e A ^] and using the Log-Sobolev inequality 
yields 

If we assume fiy(f) = 0, we find that X(0) = and so integrating with 
respect to A yields 

^(e A/ )<e A2 illl v /ll 2 H-/( 2ciV ). 
Using Chebyshev's inequality thus gives, for e > and A > 0, 

^(f>e)< e -^ e x2 W^f\\ 2 ^/( 2cN ) 
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and so optimizing with respect to A results with 

/i^(/>£)<e- ciVE2 /( 2 llll v ^ll-). 

Replacing / by — / gives the well-known concentration estimate, for any 
e>0, 

^(|/|>£)<2e- cJV£2 /(2||||V/|||| 2 00 ) 

for any continuously differentiable function / such that fiy(f) = 0. This 
estimate extends, modulo some extra technicalities, to Lipschitz functions 
and then ||||V/||||oo is replaced by the Lipschitz norm 

where x,y belong to M mAr2 and ||x|| denotes the Euclidean norm of x. Then, 
for all £ > 0, the following estimate holds: 

(9) ^(\f-^(f)\>e)<2e- N ^W\\l\ 
We set 

f P (X): = f(P)-m$„ 



= Tr(P(X)) - 4,m, 
with cf M = jv/.jv- / l A jvTr(P(A))d/^(A). Observing that 

d {Ai)kl Tr(P(A)) = (D i P(A)) lk , 

we find that on the closed set Aj^, fp is a Lipschitz (actually an infinitely 
differentiable) function of the entries of A G Hj\r(C) m with constant 

(||/p||^) 2 := sup ||VTrP(A)|| 2 

m 

= sup Y, T r( D kP(DkP)*) <N(\\P\\™) 2 , 
A6A^ k=l 

where we simply used that the set of N x N matrices is a C*-algebra. As a 
consequence, we also find that, for B £ AjJJ-, 



I/p(B) 



Tr(P(B)) - -j^l^ J t K Tr(P(A))d^(A) 



AN 

TiP 



a m / m \ 1/2 



(10) < L-Jf- J 1 A * £ Tr (^ " d$ (A) 



< 2yJmMN\\P\\f 
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and so on we can extend fp to P G C(Xi, . . . , X m )^ . 

We can also extend fp to the whole space 7i.N(C) m with the same Lips- 
chitz constant by putting 

{/ m \ l / 2 \ 

/ P (B) - \/iV||P||f £ ^ " ■ 
Vi=i / J 

Then applying (9), with 

£ p,m = I^\/( 1 (a^) C /p)I + 1 1 - mt/IA-m) -1 !!^!/ 0-a^/p)I> 

we obtain 

/#({|*V) - m£ M | > e + n A&) 



/p "^(a1)^ (1a - /p) 



> e + e ^ M nA^ 



<^(\fp-^(fp)\>e) 

A N 

< 2e - 7Vc£2 /(2(ll/p|| £ M ) 2 ) = 2e - C£ 2 /(2(||P||f ) 2 )_ 

We now use the exponential decay of the largest eigenvalue to control Ep ^ . 
By (10) and the definition of fp, note that 

/ / m \ V 2 \ 

f P (A)<2^RMN\\P\\ 1 g + ^/N\\P\\^[ K^Tr(A 2 ) + VmNM\. 

Consequently, if M,N are large enough so that //y((A^) c ) < e ~ aNM < _ ; 
by Property 2.2, 

£^ M </i^(l (A£) c(MiVm||P||f(3 + A max (A)))) + 2MrniV||P||fe- aJVM 
< M7Vm||P||^ (^5e- Q7VM + ^({A max (A) > y V M}) 



<6m( M 2 + — JiVllPl 

\ CKiV / 



M-aNM 
C e 



When P is a monomial of degree d, 

\\P\\c < dM*- 1 . 
Thus, we only need to control nip ^ ; 

< 2Ne~ aMN M d + N^(1 (a n )c A max (A) d 



+ |/^(W )c Tr(P))| 
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roo 

= 2Ne~ aMN M d + dN / y d " Vy ({Amax(A) > y V M}) dy 
Jo 

roc 

< 2Ne~ aMN M d + dN / y d-i e -aNyvM dy 

Jo 

< (2 + l)Ne~ aMN M d + dNe~ aNM V — • • • — - 

fc=i 

<iV(3M d + (i 2 )e- QjVM , 
where we assumed that d < aN. □ 



For later purposes, we have to find a control on the variance of p," . Recall 
that f{P) = N{fi N (P) - JZ N (P)). 

Lemma 2.4. For any e,rj,c> 0, there exists B, C, Mq > such that, for 
all t 6 Bjj^c , all M > Mq, for all iV£N, and all monomial P of degree less 
than eN^I^, 

${{f{P)f) < B{\\P\\^) 2 + C d N 2 e~ aMN / 2 . 
Proof. If P is a monomial of degree d, we write 

/^((IV)) 2 ) < /#(1 A * (IV)) 2 ) + m£(1(a£ )c (!V)) 2 ) = /i +h- 

(11) 

For Ji, the previous Lemma implies that 



h 



= 2 / x$({\Tr(P) - $(Tr{P))\ > x} n A£) 

/*00 

<(^ M + |m^ M |) 2 +4 / xe~ cx 2 /(2(||P||^) 2 ) dx 

Jo 



<Ce- QM ^ + i?(||P||f) 2 , 

with a constant B = - and a constant C such that (sp M + I^^m!) 2 — 
Ce -aMN f()r all d < £^2/3^ For t h e seconc i term, we take M > Mq with Mq 
as in Lemma 2.2 (exponential tail of the largest eigenvalue) to get 

h < ^[{Ki) C \ ll2 ^{{f{P)) A ) 1/2 < e-° MN / 2 ^((f(P))r /2 . 
By the Cauchy-Schwarz inequality, we obtain the control 

/^[iV) 4 ] < 2Vvf((^A Ar (p)) 4 ) < 2 4 ivVv((A J vn) 2 )- 

Now, by the noncommutative Holder's inequality (see, e.g., [21]), 

[^(PP*)] 2 < max /^(Xf) 

Ki<m 
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so that we obtain the bound 

^f{Pf]<2^N\^Jl N (Xt d ). 

l<i<m 

By (5) and (6), we obtained a uniform bound x(= 4\/ Ac" 1 ) on fi N (Xi) so 
that we have proved using (4) that 

V {Xi ) < 2 (/x c (/i (Aj )) + x ). 

We can now use the control on the moments as obtained, for instance, by 
Soshnikov (Theorem 2, page 17 in [22]) to see that there exists C(e), C(e) < 
oo for e > 0, so that 

provided d < eN 2 / 3 . As a consequence, we get that 
(12) V- N {Xf d ) < C(e) M 

for all d < eN 2 / 3 and all integer number N. Here C(s) denotes a finite 
constant depending only one, r; and c which may have changed from line 
to line. Hence, we conclude that 

h < AN 2 e- aMN ' 2 C(e) 2d . 

Plugging back this estimate into (11), we have proved that for N and M 
sufficiently large, all monomials P of degree d< eN 2 / 3 , all t G Br, tC , 

^((5 N (P)) 2 ) < B(\\P\\ h l) 2 + c 2d N 2 e~ aMN / 2 

with a finite constant C depending only on e, c and rj. □ 

3. Bound on the distance between /x and JZ N . We here bound, for all 
monomial P, 

5 N (P) = N(-p N (P)-v)(P). 



Proposition 3.1. For all c, e > 0, there exists rj > 0, C < +oo, such 
that for all integer number N , all t € B^^, and all monomial functions P of 
degree less than eN 2 / 3 , 

\s N (P)\<—. 

In particular, |(5 — 5 N )(P)\ < c ^ } almost surely. 
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Proof. The starting point is the finite dimensional Schwinger-Dyson 
equation that one gets readily by integration by parts (see [14], proof of 
Theorem 3.4) 

(13) ^(fi N [(X i + D i V)P])=^(fi N ®fi N (diP)). 
Therefore, since \x satisfies the Schwinger-Dyson equation (2) 

(14) fj,[(X i +D i V)P}=n® f J,(d i P), 

by taking the difference, we get that for all polynomial P, 

(15) 5 N (X l P) = -5 N (D l VP) + 5%7Z JV (d i P) + »®5 N (diP) +r(N,P), 
with 

r(N,P) :=N- 1 ^(f ®f(diP)). 

If we take P a monomial of degree d < eiV 2/3 and assume M > Mo, then we 
see, by using Lemma 2.4, 

wn,p)\<± e ^(\f(Pi)\ 2 ) 1/2 ^m N (P2)\ 2 ) 1/2 

iv p = p lxi p 2 

C d-1 



< ^Y(Bl 2 M 2 «-V + c l N 2 e- aMN / 2 ) 1/2 
- N z^y 



1=0 



We set 



x ( B (d - I - ifM 2 ^- 1 ^ + ^(^1-1)^6-^/2)1/2 

< ^d{B{d-l) 2 M 2{ - d - 2 ) +C ^N 2 e- aMN / 2 ) :=r(N,d,M). 



A% = max \S N (P)\. 

P monomial of degree d 



Observe that by (12), for any monomial of degree d less than eiV 2 / 3 , \fi N (P)\ < 
C(e) d , \n(P)\ <C$< C(e) d . It allows us to obtain the rough bound < 
2NC(e) d if d < eN 2 / 3 . By (15), writing DiV = E^Agy, we get that, for 

d<eN 2 / 3 , 



d-l 



A^ +1 < max £ \ tj | A» +deg{Diqj) + 2 £ C{e)^ A? + r(N, d, M). 



3=1 1=0 

We next define, for k < 1, 

A N ( K ,e):= J2 ^k- 

k=l 
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We obtain, if D is the maximal degree of V, 

A n (k,e) < [C' K ~ D \t\ + 2(1 - C(e)Ky l K 2 ]A N ( K ,e) 

(16) 

£ N 2 / 3 +D eN 2 / 3 

+ C\t\ £ K k ~ D ^k+ E * h+1 r(N,k,M), 

k= £ N 2 / 3 + l fc=l 

where we choose k small enough so that C(e)k < 1. Moreover, since D is 
finite, using the bound on A^f , we get 

£ N 2 / 3 +D 

£ K h - D ^ < 2DN(kC(s)) £n2/3 k- d . 

k=eN 2 /' 3 +l 

Since kC(e) < 1, as N goes to infinity, this term is negligible with respect 
to iV -1 for all e > 0. The following estimate holds: 

eN 2 / 3 

]T K k r(N,k,M) 
k=l 

n eN 2 /' 3 n n 

< - £ kn k {B{k - l) 2 M 2 ( fc ~ 2 ) + c^N 2 e~ aNM l 2 ) < — 
k=i 

if k is small enough so that M 2 k < 1 and Ck < 1. We observed here that 
j\j2 & -aNM /2 j g un iformly bounded independently of N £ N. Now, if |t| is 
small, we can choose k so that 

C := 1 - [C'K- D |t| + 2(1 - C(e)k)- 1 k 2 \ > 0. 

Plugging these controls into (16) shows that for all e > 0, and for k > small 
enough, there exists a finite constant C(k, e) so that 

A JV (K,e)<C r (K,e)JV- 1 

and so for all monomial P of degree 

|tf Ar (P)|<C f (K,e)K- d i\r 1 . □ 

To get the precise evaluation of (P) , we shall first obtain a central limit 
theorem under fly which in turn will allow us to estimate 

lim Nr(N,P). 

N^oo 
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4. Central limit theorem. We shall here prove that 

S N (P)=N(/1 N -^(P) 

satisfies a central limit theorem for all polynomial P. By Proposition 3.1, it 

is equivalent to prove a central limit theorem for 5 (P) , P G C{X\ , . . . , X m ) . 
We start by giving a weak form of a central limit theorem for Stieljes-like 
functions. We then extend the result to polynomial functions in the image of 
some differential operator. We finally generalize our result to any polynomial 
functions. 

For the rest of the paper, we will always assume the following hypothesis 
(H). 

(H): Let c be a positive real number. The parameter t is in B VjC with r] 
sufficiently small such that we have the convergence to the solution of (2) as 
well as the control given by Lemma 2.1 (Compact support) and Proposition 
3.1. 

Note that (H) implies also that the control of Lemma 2.1 (Compact sup- 
port) is uniform, and that we can apply Lemma 2.2 (Exponential tail of the 
largest eigenvalue) and Lemma 2.3 (Concentration inequality) with uniform 
constants. 



4.1. Central limit theorem for Stieljes test functions. One of the issues 
that one needs to address when working with polynomials is that they are 
not uniformly bounded. For that reason, we will prefer to work in this section 
with the complex vector space C^(C) generated by the Stieljes functionals 

(17) ST m (C) = l J] f-Zi-E^Xfc) ; ^GC\R,a?€R,p€Nl, 
ll<i<p\ k=l ) J 

where is the noncommutative product. We can also equip ST m (C) with 
an involution 

/ -» / m \ _1 \ * -+ / m \ -1 

n u-e^x, = n fc-E«r fcx • 

\l<k<p\ i=l ) ) l<k<p\ i=l J 

We denote C^(C) sa the set of self-adjoint elements of C^(C). The derivation 
is defined by the Leibniz rule and 

(m \ ~ 1 / m \ / m 

We recall notation; first jj is the operator 




(P®Q)$h = PhQ and (P <g> Q ® R)${g,h) = PgQhR 
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so that, for a monomial q, 

di° djq#(hi,hj) = ^ qohiqihjq 2 + ^ qohjqihiq 2 . 

q=qoXiqiXjq 2 q=qoXjq 1 X i q 2 

Lemma 4.1. Assume (H) and Zei h±,...,h m be in C^(C) sa . Then the 
random variable 

m 

Y N (h u ...,h m ) = Nj2i^ N ® - + A^fc]} 

fc=l 

converges in law toward a real centered Gaussian variable with covariance 

m m 

C{h l ,...,h m )= (v®v[d k h l xdih k }+ii(diod k V$(h k ,hi))) + J2Kht)- 

k,l=t k=l 

PROOF. Define W = \ Y,i X i+V- Notice that Y N (hi, . . . ,h m ) is real val- 
ued because the h' k s and W are self adjoint. The proof follows from the usual 
change of variable trick. We take hi, ... , h m in C™(C) sa , A £ R and perform 
a change of variable Ai — > Bi = F(A)i = A{ + -4/jj(A) in Zy . Note that since 
the hi are C°° and uniformly bounded, this defines a bijection on 7Yjv(C) m 
for iV big enough. We shall compute the Jacobian of this change of variables 
up to its second order correction. The Jacobian J may be seen as a matrix 
(Ji. j)i<i,j<m where the Jij are in Z^(7Yjy(C)) the set of endomorphisms of 
TCn(C), and we can write J = I + -4 J with 

7 M -:W 7V (C)^H A r(C), 
X — ► dihj#X. 

Now, for 1 <i,j <m, X — > dihj^X is bounded for the operator norm uni- 
formly in N [since hj S C s t(C),dihj S C s t(C) <8>C s t(C) is uniformly bounded] 

v 



so that, for sufficiently large N, the operator norm of -^J is less than 1. 



From this, we deduce 

det / + - 



det J| 



Observe that as J is a matrix of size m N and of uniformly bounded norm, 

the kth term - — ^ Tr(J ) is bounded by N l_ 2 . Hence, only the two 

first terms in the expansion will contribute to the order 1 and the sum sn 
of the other terms will be of order -k. To compute the two first terms in 
the expansion, we only have to remark that if (ft is an endomorphism of 
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H N (C) is of the form 4>(X) = ^ AjXBj, with N x N matrices Ai,B h then 
Ti<f> = J2i TrA; TrBi [this can be checked by decomposing on the canonical 
basis of Hjv(C)]. Now, 

1% dihi 2 m 2 h i3 t(- ■ ■ (VA-PO • • •))• 

Thus, we get 

Tr(J) = ^TrJ ii= ]T Tr ® Tr(c^) 

i l<i<m 

and 

Tr (^ 2 )=E Tr ( 7 ^)= E Tr<g,Tr(c^6^)- 
We now make the change of variable Ai — > + -^h(A) to find that 

z JV = | e - J VTr(y(A)) d/iJ V (A) 

e -JVTV(W(A i +(A/JV)/ li (A))-W( J 4 i )) e (A/iV)X; i Tr®Tr(9 j ? l i) 



x e (-A 2 /(27V 2 ))^ ij Ti-®T t (a i h. J 9 J ^) eSivd/l 7V^ A ) 

with stv of order The first term can be expanded into 

w U + ^) - w(A) = 1 £ am + ^,o a,w#(k, M + ||, 

where -Rat is a polynomial in the /ij's and in the AYs, of degree less than the 
degree of V minus two in the later. To sum up, the following equality holds: 

e XY N (h 1 ,...,h m )-(\ 2 /2)C N (h 1 ,...,h m )+(l/N)(fi N (R N )+Ns N ) = 1 



with 



C N (ht, ...,h m ):= fi N (j^ d < ° d jW#(h h hj^j + A* ® A* (j2 9 i h J d i h ^j ■ 
We can decompose the previous expectation in two terms E\ and E2 with 

E 1 =Hy[l A N e XY N(hi,.-.,h m )-(\y2)C N (h 1 ,...,h m )+(l/N)(il N (R N )+Ns N )^ 

and 

E 2 = j^N^^ ^ e XY N (hi r ..,h m )-(X 2 /2)C N (hi r ..Am)+{l/N)^ 
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We first consider E\. On A^j = {A : max,,(A^ ax ( J 4j)) < M} the polynomial 
Rn is uniformly bounded and so P, n (Rn) + Nsn is of order one, bounded 
by a constant which goes uniformly to when N goes to infinity. We 
next show that we can replace Cjv(/ii, . . . , h m ) by its limit C(h\, . . . , h m ) in 
the exponential in E\. An intermediate step is to replace it by 



C N (hi,...,h m ) =^ 



\ i,j / \ i,j / 



In fact, by Lemma 2.3, fi N {P) converges toward its expectation ~p N (P) under 
fiy (l^iv •) except on sets with probability of order e~ N once evaluated at 
any products of the h^s and the X^s (because the Lipschitz constant of 
finite products of h^s and AVs are bounded on and the error terms 
e PM ano - m PM can b e bounded as we did for polynomials). Hence, we can 
find a constant C(M, c) > such that for iV large enough, 

l#({\C N {h u . ..,hm)- C N (hx,. ..,h m )\> 2e} n A&) 



< 2e 



-C(M,c)N 2 (s N ) 2 



with en = £ — £pM ~~ m PM ~ £ - Moreover, JI N (P) converges to n(P) for 
any polynomial function P (see [14], Theorems 3.1 and 3.4). Since by the 
Weierstrass theorem the hi 's can be uniformly approximated by polynomials 
on Ajy-, uniformly in N, we also know that C]y(hi, . . . , h m ) converges to 
C(h\, . . . , h m ). Consequently, we obtain for some positive constant C(M, c), 
iV large enough 

f$({\C N (hi,. ..,h m )- C(h u . ..,h m )\>e}n A&) 

< 2e- c{M ' c)N ^ 6N)2 . 

Finally, Yjv(/ii, . . . , h m ) is at most of order N and Cjv(/ii, . . . , h m ) of order 
one. Hence, the exponential in E\ is at most of order e CN for some finite 
constant C. Therefore, if we let 



E[:=^[l K e 



N \1\ „ p XY N (h 1 ,...,h m )-(\y2)C(h 1 ,...,h m h 



we deduce that 



W U N \1 Ne XY N(h 1 ,...,h m )-(X 2 /2)(C N (h 1 ,...,h m )-C(h u ...,h ri 

log — <- i— ^ V M 

E l 



< 



loge 



^[l A *e* y *(fci,-,ftm)] 



< \l og (e^/ 2 ^+2e CN e- c ^ N2£ N)\+A N . 
Letting first N going to infinity and then s going to zero yields 
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Note that this estimate is valid for any M large enough so that Lemma 
2.3 holds. 

Our goal is now to show that, for M sufficiently large, E2 vanishes when 
N goes to infinity. It would be an easy task if the term in the exponential 
were bounded, but it may in fact be large due to some derivatives of V 
appear so that there are polynomials term in the exponential. The idea to 
pass this difficulty is to make the reverse change of variables. For N bigger 
than the norm of the fej's, and with Bi = A4 + -^/ij(A), 

&2 — H>V [ JL {A:A JV (A)>M} e J 

= ^(B : A^ ax (A) > M) <(4(B: A^(B) > M - 1). 

This last quantity goes exponentially fast to for M sufficiently large by 
Lemma 2.2 (exponential tail of the largest eigenvalue). 
Hence, we arrive, for M large enough, at 

(18) lim / t ANe ^N(h 1 ,...,h m ) dfl N = e (xy2)C(h 1 ,...,h m )_ 

Since jiy (A^ f ) goes to one as N goes to infinity, we find that Ypj{h\, . . . , /i- m ) 
converges in law under fiy (A^-) _1 //y (• n A^) toward a centered Gaussian 
variable with covariance C(h\, . . . , h m ), for any M large enough. For the 
same reason, we conclude that the same convergence holds under fiy . □ 

4.2. Central limit theorem for some polynomial functions. We now ex- 
tend Lemma 4.1 to polynomial test functions. 

Lemma 4.2. Assume (H) and let P\,...,P m be in C(Xi, . . . ,X m ) sa . 
Then, the variable 

m 

Y N (P 1 ,...,P m ) = Nj2 \A N ® » N {dkPk) ~ pi N [{X k + D k V)P k }] 

k=l 

converges in law toward a real centered Gaussian variable with covariance 

m m 

c(p 1 ,...,p m )= J2 (/i^[^x9,p fc ]+/i(a l o9 fc ^(p^i)))+E/ i ft 2 )' 

k,l=l k=l 

Proof. Let Pi,...,P m be self-adjoint polynomials and h\,...,h^ n be 
Stieljes functionals which approximate Pi, ... , P m such as 

/lf(X) = jR '(l^F'---'T^Al)- 

Since E[Y N ] =0 by (13), 

A AT 

Y N (P 1 ,...,P m )=5 (K N (P 1 ,...,P m )), 



ASYMPTOTICS FOR MATRIX MODELS 



25 



with 

m 

K N (P U ...,P m ) = J2(P' N ® HdkPk) ~ (X k + D k V)P k ) 
k=l 

-TV 

and, similarly, Y^{h\^ . . . ,h £ m ) = 5 {K^{h\ , . . . , h £ m )). It is not hard to see 
that 

\\K N (hl,. ..,h £ J- K N (P U P m )\\% < eC{M) 

for some finite constant C(M) which only depends on M. Hence, we de- 
duce by Lemma 2.3 (Concentration inequality) that there exists mp £ M and 
Ep £ M going to zero as N goes to infinity (note here that the control on 
m Pe m an d £ p £ m follows exactly the same line as for monomials) such that 

^(\5 N (K N (h\, . . . , h%) - K N (P 1: ...,P k ))- m% M \ > 5 + e N P ^ M ) 

< e -aMJV + e -«5 2 /(2 C£ 2 C(A/) 2 ) 

and so for any bounded continuous function / : R — > R, if v a i is the centered 
Gaussian law of covariance a 2 , 

lim ^(f(5 N (K N (P 1 ,...,P k )))) = ]im lim l 4{f{f{K N {h\,...,h%)))) 

TV— >oo e— >0TV— too 



limu Cih[ ,..., h e m) (f) = u C ( Pl ,..., Pm) (f), 



where we used in the second line Lemma 4.2 and in the last line Lemma 2.1 
(Compact support) to obtain the convergence of C(h\, ■ . ■ , h^) to C(P\, . . . , 
Pm). □ 

Yjv depends on Nfi N <S> ft , in which clearly one of the empirical distribu- 
tion fi N shall converge to its deterministic limit. This is the content of the 
next lemma. 

Lemma 4.3. Assume (H) and let Pi, . . . ,P m be self-adjoint polynomial 
functions. Then, the variable 

(mm \ 
]T (x k + D k v)p k -^2(h®i + i®h) (d k p k ) 
k=l k=l ) 

converges in law toward a centered Gaussian variable with covariance 

m m 

C(Pi, . . . , P m ) = (M® l*[dkPi x diPk] + Kdi o d k V^P k ,Pi))) + £ MP?)- 

k,l=l k=l 
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Proof. The only point is to notice that using (2), 

m 

Y N (P U . . . ,P m ) = Y,0 N ® M + $ N )(dkPk) ~ S N ((X k + D k V)P k ) + r NtP 
k=i 

with r^p = A^ 1 X^feLi " ® v {dkPk) of order A^" 1 with probability going 
to 1 by Lemma 2.3 (Concentration inequality) and Property 3.1. Thus, 

y*(Pi,...,p m ) 

= ^ + DkV ^ Pk + ® M + ® I){dkPk))^j + r N ,p 

= -Z N (P 1 ,...,P m ) + o(jj 
This, with the previous lemma, proves the claim. □ 

4.3. Central limit theorem for all polynomial functions. In the previous 
part we have obtained CLT's only for the family of random variables S N (Q) 
with Q in the following subset T of C(Xi, . . . ,X m ): 

{mm \ 
Y,(Xk + D k V)P k + ii)(dkPk),Vi, Pi self-adjoint . 

ib=i k=i ) 

In this section we wish to extend it to 5 N (Q) for any self-adjoint polynomial 
function Q, that is, to prove Theorem 1.3. We have to show a form of density 
oi^mC(X 1 ,...,X m ). 

The strategy is to see T as the image of an operator that we will invert. 
The first operator that comes to mind is 

m m 
: (Pi, . . • , P fc ) -> ^(X fc + L^P* - ]T> ® / + 1 <8> M)(3fcPfc) 
fc=l fe=l 

as we immediately have T = ^(C(Xx, . . . ,X m ) sa , . . . , C(Xi, . . . ,X m ) sa ). 

In order to obtain an operator from C(Xi, . . . , X m ) sa to C(X±, . . . , X m ) sa , 
we will prefer to apply $ to = -DfcP for all k and for a given p as we 
shall see later, ^{D\P, . . . , D m P) is closely related with the projection on 
functions of the type TrP of the operator on the entries A — VA r Tr(Vl^).V 
which is symmetric in L 2 ([iy). The resulting operator is a differential oper- 
ator and, hence, it would be hard to prove that it is continuous on a fixed 
space of functions. To avoid this issue and make the argument more readable 
we have first to divide each monomial of P by its degree. 

Then, we define a linear map S on C{X\, . . . , X m ) such that, for all mono- 
mials q of degree greater or equal to 1, 
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Moreover, S(g) = if degq = 0. For later use, we set Co(-X"i, • • • ,X m ) to be 
the subset of polynomials P of C{Xi, . . . , X m ) sa such that P(0, . . . , 0) = 0. 
We let n be the projection from C(Xi, . . . , X m ) sa onto Co(-X"i, ■ • ■ , X m ) [i.e., 
n(i- > ) = P — P(0, . . . , 0)]. We now define some operators on Co(Xi, . . . , X m ) 
that is, from Co , ■ • ■ , X m ) into Co (X\ , • ■ • , X m ) , 

E 1 :P^n(^2d k ZP$D k V ) J, 

E 2 :P^U ® J + 7® /J,)(d k D k XP)j . 

We denote E = Id — S 2 and S = S + Si, where I is the identity on 
Cq(Xi, . . . , X rn ) . Note that the images Sj's and S are indeed included in 
C(-Xi, . . . , X m ) sa since V is assumed self-adjoint. With this notation, Lemma 
4.3, once applied to Pi = DiT,P, 1 < i < m, reads as follows: 

Proposition 4.4. For all P in Co(X\, . . . ,X m ), 5 N (EP) converges in 
law to a centered Gaussian variable with covariance 

C(P) :=C(D 1 EP,... ) D m EP). 



Proof. We have for all tracial state r, T(d k P$V) = r(D k PV) and if P 
is in Co(A^i, . . . ,X m ) [i.e., P(0, ... ,0) = 0], we have the identity 

P = Y,d k ^P%X k . 

k 

Then, as 8 is tracial and null on constant terms (so that the projection II 
can be removed in the definition of S), for all polynomial P, 

(m m 
P + Y, d k XP$D k V - + ii){d k D k -ZP) 

k=l k=l 

(m m 
(X k + D k V)D k EP -]T(/i®H-I® fi)(d k D k XP) 
k=l k=l 

= Z N (D 1 ZP,...,D m ZP). 
We then use the Lemma 4.3 to conclude. □ 



To generalize the central limit theorem to all polynomial functions, we 
need to show that the image of S is dense and to control approximations. If 
P is a polynomial and q a nonconstant monomial, we will denote X q (P) the 
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coefficient of q in the decomposition of P in monomials. We can then define 
a norm || • \\a on Co(^i, • • • ,X m ) for A > 1 by 

||P|U= Yl \K(P)\A desq . 

deg q^O 

In the formula above, the sum is taken on all nonconstant monomials. 
We also define the operator norm given, for T from Co{Xi, . . . , X m ) to 
C Q (X 1 ,...,X m ), by 

|||T||U= sup \\T(P)\\ A . 
\\p\\a=i 

Finally, let Co(^i, • • • ,X m )A be the completion of Co(Xi, • • • ,X m ) for || • \\a- 
We say that T is continuous on Co(Xi, • • • ,X m )A if \\\T\\\a is finite. We shall 
prove that H is continuous on Cq(Xi, . . . , X rn )A with continuous inverse when 
t is small. 

Lemma 4.5. With the previous notation: 

1. The operator Ho is invertible on <Cq(Xi, • • • ,X m ). 

2. There exists Aq > such that, for all A > Aq, the operators H2, Ho and 
Hq 1 are continuous on Cq(X\, . . . , X m ) a and their norms are uniformly 
bounded for t in B v . 

3. For all e,A > 0, there exists n £ > such for |t| < rj £ , Hi is continuous on 
C {Xi, . . .,X m ) A and \\z>i\\a < e. 

4. For all A > Aq, there exists r/ > such that for t £ B v , H is continuous, 
invertible with a continuous inverse on Cq{X\, . . . ,X m ) a- Besides, the 
norms o/H and H" 1 are uniformly bounded for t in B^. 

5. There exists C > such that, for all A> C, C is continuous from Co(^i , • • • , 
X m )A into 

Proof. 1. We can write 

Ho = I — H2. 

Observe that since H2 reduces the degree of a polynomial by at least 2, 

p-£(h 2 hp) 

n>0 

is well defined on Co(A^i, . . . , X m ) as the sum is finite for any polynomial P. 
This clearly gives an inverse for Ho- 

2. First remark that a linear operator T has a norm less than C with respect 
to || • \\a if and only if for all nonconstant monomial q, 

\\T(q)\\ A <CA dc ^. 
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Recall that \x is uniformly compactly supported [see Lemma 2.1 (Compact 
support)] and let Co < +00 be such that \p(q)\ < C egq for all monomial 
q. Take a monomial q = ■ ■ ■ X;, and assume that A > 2Cq, 



nlj2(I®v)dkD k Xq 



<P 1 E ll r iAt(r 2 )|U 

k,q=qiX k q 2 , 
q2qi=riX k r2 

p-lp-2 

< p- 1 y, ^ dcgri c dcgr2 = - E E a 1 ci~ 1 - 2 

k,q=q^X k q 2 , ^ n=0 1=0 



qiq\=r\X k ri 
P-2 / r~i \ p-2-l 

<A p - 2 Y,(— 



<2A- 2 \\q\\ A , 



where in the second line, we observed that once deg(gi) is fixed, q 2 qi 
is uniquely determined and then r\,r2 are uniquely determined by the 
choice of I the degree of t\. Thus, the factor - is compensated by the 
number of possible decompositions of q, that is, the choice of n the degree 
of q\. If A > 2, P — > II(^)j fc (7 <8> fJ^df-D^HP) is continuous of norm strictly 
less than |. And a similar calculus for n(J^j fe (/i (8) 1)8^0^) shows that 
H2 is continuous of norm strictly less than 1. It follows immediately that 
Hq is continuous. Recall now that 



"0 



n>0 



As H2 is of norm strictly less than 1, Sq 1 is immediately continuous. 
3. Let q = X^ ■ ■ ■ Xi p be a monomial and let D be the degree of V and 
B{< Dn) the sum of the maximum number of monomials in D k V: 

< - E \\qiD k Vq 2 \\A 



P 



k,q=qiX k q 2 



< 



1 

P 



E 



\t\BA p 



-l+D-l 



k,q=q\X k q 2 

= \t\BA D - 2 \\q\\ A . 
It is now sufficient to take r\ e < (BA D ~ 2 )~' 



s. 



4. We choose 77 < (BA D 2 ) 1 |||H 1 ||| A 1 so that when |t| < 77, 
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By continuity, we can extend Ho, Hi, H2, E and Hq 1 on the space Co(-Xi, . . . , 
X m )A- The operator 

n>0 

is well defined and continuous. And this is clearly an inverse of 

H = H + Hi = H (7 + H 1 Hi). 

5. We finally prove that C is continuous from Cq(Xi, . . . , X m )A into R where 
we recall that we assumed A > Cq. Let us consider the first term 

m 

d(P) := A*®M(%Dl£P x dtD k EP). 

k,l=l 

Then we obtain, as in the second point of this proof, 

9291 =n X; r 2 ,q' 2 q[ =r[ X k r' % 

<4^|A,(P)||A (? ,( J P)|deggdeg (? 'C deg9+des,? '- 4 
9,9' 

<4fsup^- 2 ^Vlli- 

\£>0 / 
We next turn to show that 

m 

c 2 (P) := X! K^°<W(£>feEP,A£P)) 
fc,/=i 

is also continuous for || • \\a- In fact, noting that we may assume V € 
Co(-Xi, . . . , X m ) without changing C2, 

|C 2 (P)|< J2 \X P (V)\ 

p,q,q',k,l 

I A g (P) ||A 9 /(P)| c dcgp+dcg 9+deg 9 '" 4 



x 



y 

, v v deggdegg' 

q,q' ,p=piX k p 2 Xip 3 
q=qiX k q2,q'=q[X k q' 2 

< n\t\D 2 Y / \K( p )\\^q'( p )\Co +dcsq+dcsq '~ A 



9>9 

ii-nu- 



<n|t|P 2 G?- 4 ||P" 2 



The continuity of the last term C%{P) = YaLi ^((Dj^P) 2 ) i s obtained 
similarly. □ 
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We can compare the norm || • \\a to a more intuitive norm, namely, || • \\^ 
denned in (7). 

We will say that a semi-norm M is weaker than a semi-norm N 1 if and 
only if there exists C < +00 such that, for all P in Co(-Xi, • • • , X m ), 

Af(P) < CN'{P). 

Lemma 4.6. For A > M , the semi-norm \\ ■ \\^ restricted to the space 
Cq(X±, . . . , X m ) is weaker than the norm \\ ■ \\a- 

Proof. For all P in Co(-Xi, • ■ • ,X m ), the following inequalities hold: 
\\P\\% < E MP)\U\\c < E I V)l deggM de ^ < (supi(^)') I^IU b 

To take into account the previous results, we define a new hypothesis (H') 
stronger than (H). 

(H'): (H) is satisfied, A — 1 > max(Ao,Mo,C) for the Mo which appear 
in Lemma 2.2 (Exponential tail of the largest eigenvalue) and the C which 
appear in Proposition 3.1. Besides, |t| < 77 with 7/ as in the fourth point of 
Lemma 4.5 in order that 3 and H _1 are continuous on Co(X\, . . . , X m )A and 
Co(-Xi, . . . ,X m )A-i, and that C is also continuous for these norms. 

The two main additional consequences of this hypothesis are the conti- 
nuity of S for || • The strange condition about the continuity of H on 
Co(Xi, . . . ,X m )A-i is here for a technical reason which will appear only in 
the last section on the interpretation of the first order correction to the free 
energy. 

While (H') is full of conditions, the only important hypothesis is the 
convexity of V. Given such a V, we can always find constants A and rj which 
satisfy the hypothesis. The only restriction will be then that t is sufficiently 
small. 

We can now prove the general central limit theorem which is up to the 
identification of the covariance equivalent to Theorem 1.3. 

Theorem 4.7. Assume (H'). For all P in C(X U . . . , X m ) sa , 5 N (P) con- 
verges in law to a centered Gaussian variable 7p with covariance 

a\P) :=C{~- l Tl(P)) = CiPxYEr^P), . . .,D m TSr x U{P)). 

If P € C(Xx, ■ ■ ■ ,X m ), 5 N (P) converges to the complex centered Gaussian 
variable "f(p + P*y2 + i^frp-p*)/2i [the covariance of 7(p+p*)/2 and 7(p_p*)/2i 
being given by a 2 ((P + P*)/2, (P — P*)/2i), where <r 2 (-, •) is the bilinear form 
associated to the quadratic form a 2 ]. 
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Proof. As 5 n (P) does not depend on constant terms, we can directly 
take P = n(P) in Co(^i, • • • ,X m ). Now, by part 4 of Lemma 4.5, we can 
find an element Q of Co(Xi, . . . , X m ) A such that EQ = P. But the space 
Co(Xi, . . . , X m ) is dense in Co(X\, . . . , X m ) A by construction. Thus, there 
exists a sequence Q n in Co(Xi, ■ • ■ , X m ) such that 

lim \\Q-Qu\\a = 0. 

n—too 

Let us define R n = P — EQ n in Co(-Xi, • • • , X m ). 

Now according to Property 4.4 for all n, 6 (HQ„) converges in law to a 
Gaussian variable 7„ of variance C(Q n ) with 

C(Q n )=C(D 1 ZQ n ,..-,D m ZQ n ). 

As C is continuous by part 4 of Lemma 4.5, it can be extended to the 
space C {X 1 ,..., X m ) A and a 2 (P) = C(H~ 1 P) = C(Q) = lim n C(Q n ) is well 
defined. Hence, 7 n converges weakly toward 700, the centered Gaussian law 
with covariance C(Q), when n goes to +00. The last step is to prove the 
convergence in law of 5 N (P) to 700- We will use the Dudley distance. For 
/ : M — ► R, we define \f\c = \\f\\c + ||/||oo- The Dudley distance between two 
measures on R is 

V(n,u)= sup \v(f)-v(f)\. 

l/k<i 

The topology induced by the Dudley metric is the topology of the conver- 
gence in law. Below, as a parameter of V, we denote in short 5 N (P) for the 
law of 5 N (P). We make the following decomposition: 

V(S N (P), loo )<V(S N (P),S N (EQ n )) 

(19) 

+ V(S N (EQ n ), 7n ) + £>( 7n , 7oo ). 

By the above remarks, T>(5 N (EQ n ), j n ) goes to when N goes to +00 and 
P(7 n ,7 oc ) goes to when n goes to +00. We now use the bound on the 
Dudley distance: 

V(5 N (P),5 N (EQ n )) < E[\5 N (P) - 6 N (EQ n )\ A 1] = E[\5 N (R n )\ A 1]. 

We control the last term by Lemmas 2.3 (Concentration inequality) and 2.2 
(Exponential tail of the largest eigenvalue) so that, for M > Mo, 

E[\5 N (R n )\ A 1] < e -«™ + 2^||i? n ,||f + el M + \m N RnM \. 

But we deduce from Lemma 4.6 that since we chose M < A, there exists a 
finite constant C such that 

IKH^ < C\\R n \\ A = C\\E(Q - Q n )\\ A < C|||S||m|Q - Q n \\ A 
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and so H-Rnll^ g° es to zero as n goes to infinity. And since ||-Rn||£ * s finite, 
e R n M g° es to zero. Similarly, using the bound of Lemma 2.3 on rrip M for 
P monomial, we find that 

|ml, M l <^E \X q (R n )\deg(q)(3M d ^ + deg(q) 2 )e- aMN 
q 

< N sup{£(3M e + £ 2 )A- e )\\R n \\ A e- aMN 

goes to zero as ./V goes to infinity. Thus, E[\5 N (R n )\ A 1] goes to zero as n 
and N go to infinity. Putting things together, we obtain if we let first iV 
going to +oo and then n, the desired convergence lim.]yV(5 N (P), 700) = 0. 
□ 



Note that the convergence in law in Theorem 4.7 can be generalized to a 
convergence in moments; 

Corollary 4.8. Assume (H 7 ). Let P be a self-adjoint polynomial, then 
S N (P) converges in moments to a real centered Gaussian variable with vari- 
ance a 2 (P), that is, for all k in N, 

hm / (S N P) k d$ = , \ >/ / zV* 2 /^)) dx. 

Proof. Indeed, once again we decompose f(5 N P) k d^ into + E? 
with 

< = / ± AN J N Pfd$E$ = J t [K)c (8 N P) k d^, 

with M > Mq. For E±, we notice that the law of 5 N P has a sub-Gaussian tail 
according to Lemma 2.3 (Concentration inequality). Therefore, we can re- 
place x k by a bounded continuous function, producing an error independent 
of N . Applying Theorem 4.7 then shows that 

hm / t AN (5 N P) h d$ = ; * ; / ^e-* 2 /(2- 2 GP)) dx . 
For the second term, we use the trivial bound 

\E?\<N k J l (A ^ )c (|A max (A)| + H(P)) fc ^ 

< kN k f (A + \n\ (P)) k ~ 1 e- aXN dX, 

J\>M 

which goes to zero as N goes to infinity for all finite k. □ 
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Another generalization of Theorem 4.7 is to extend the set of test functions 
from polynomials to the completion of Cq(X%, . . . ,X m ) for the Lipschitz 
semi-norm || • ||^ . We shall assume that M is strictly greater than C, the 
constant which bounds uniformly the radius of the support of /i according 
to Lemma 2.1 (Compact support), and also greater than M , the constant 
which appears in Lemma 2.2 (Exponential tail of the largest eigenvalue) 
in order to have A max (A) less than M with high probability. We denote 
Co(-X"i, • • • ,X m )£ the completion of Co(-Xi, ■ • ■ ,X m ) for that norm. 

Let us first extend some of the previous quantities to this setting. Recall 



that, for all TV G N, 
A 



N\\P\\% 



is always bigger than ||TrP 



\ M 



so that if 



(A) < M, TrP(A) is well defined. This allows us to define, for P in 
Co{Xi,. . .,X m )%, fi N {P) = ^TrP(A) on Aj£. We can also extend p to this 
context by the following: 



Lemma 4.9. Let P £ Co(-Xi, • • • ,X m ). Then, with Co as in Lemma 2.1 
(Compact support), 

MP)\<V^C \\P\\Z°. 
Proof. Let us consider the following norm on C(X±, . . . ,X m ): 
||P|| M :=limsup(M(PPT)) 1/(2n) . 

n 

The completion and separation of C{X\, . . . ,X m ) for this norm is then a 
C*-algebra (see, e.g., the Gelfand-Neimark-Segal construction). As \i is 
compactly supported, the norm of the Xj's are bounded by Co- Besides, 
for all P, 



Therefore, we can write 

Mp)I = Kp(x))-m(p(o))| 



/ J2(Dk p )( s ^) x kds 



< 



Jj^D k P(sX)X k 



<k=l 



ds 



^ l (£M^(*X)L> fe P( S X)*)J f$>(*g)J ds 



1/2 



I Co 



<C sup Ell- D fc- P ( a; i)---> x m)||^ =VrnC \\P\ lc . 

„4C*-algebra \fc=l / 
a;i=a:*||^i||<Co Q 
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Thus, \x extends to C(X X , X m )f . It is a natural question to study the 
behavior of 

5 N (P) ^N^iP)-^))^ 
for P in Co{X u . . . , X m )]f, the completion of C {Xx, X m ) for || • ||f . 

Corollary 4.10. Assume (H') and let M be bigger than Cq and M : 

1. a 2 is continuous for || • and so extends to Cq(Xi, . . . ,X m }^ . 

2. For all P in Cq{X\, . . . ,X m )^ , 5 N (P) converges in law to a Gaussian 
variable with variance a 2 (P). 

Proof. We take a sequence of polynomials S n which converges to P for 
the norm || • \\^ . Let R n = P — S n be the rest. For all n, 5 N (S n ) converges 
to a centered Gaussian variable 7 n of variance a 2 (S n ). 

Let us show that a 2 is continuous for || • . Let P be a polynomial, and 
M sufficiently large, 

a 2 {P) = lunE[f(P) 2 } =limE[t K 8 N (Pf]. 

The first equality comes from the previous corollary about the convergence 
in moments, as well as Lemma 3.1, which allows to recenter with respect 
to the mean rather than the limit, and the second equality comes from 
Lemma 2.2 (Exponential tail of the largest eigenvalue). Now by Lemma 2.3 
(Concentration inequality), as ||-P||/f controls the Lipschitz norm of iTr(P), 

lim M £[WlV)) 2 ] 

N— >oo M 

= 2 lim / erfKA&nte (P)\>s})de 

< r2ee^ 2 l^c) 2 ) d e = -{\\P\\^) 2 , 
Jo c 

where we used that mj^- P and ej^-p of Lemma 2.3 go to zero as N goes to 
infinity since P is a polynomial. Thus, the quadratic form a 2 is continuous for 
|| • ||£ f and can be extended on Co(-X"i, • • • ,X m )^ . This implies that a 2 (S n ) 
converges to a 2 (P). The rest of the proof is exactly as that of Theorem 4.7 
and we omit it. □ 

Note that by Lemma 4.5 the norm || • \\a is stronger than the norm || • 
so that we can use this corollary to extend out the central limit theorem on 
Cq(X±, . . . ,X m )A and, by continuity of a 2 , on this space the formula 

a 2 (P) := C(H" 1 P) = CiD^E^P, . . . , D^EE^P) 

remains valid. 
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5. Identification of the variance. 

5.1. Exact formula. We shall provide here a more tractable formula for 
the variance o~ 2 (P) of the limiting Gaussian distribution found in Theorem 
4.7. Note that for all polynomials P, Q, 5 N (P + Q) converges to Jp+q- 
Thus, {7p|Pg C(Xi,...,X m ) sa } = {j P \P£ C {Xi,...,X m )} has a natural 
structure of Gaussian space. In this space all elements are centered and the 
covariance function is given, for P,Q G Co(-Xi, • • • , X m ) by 

a\P, Q) =C(Z~ 1 P, E^Q) = C(D£~^P, DES^Q), 

where D is the cyclic gradient defined by DP = (AP . . . , D m P) and 

C(P\, . . . , P m , Qi, . . . , Qrn) 

m m 

= J2 (A* ® v[d k Pl x diQk\ + »(Pl o d k V$(P k , Q,))) + £ At(P fc Q fc ). 
fc,«=i fc=i 

We now give a more explicit formula for a 2 (P,Q). We therefore need to 
study C and the commutation relations of the cyclic gradient and E. 
Let us define the following operators on C(Xi, . . . , X m ): 

m m 

Si :P — £ dfcPttlW H 2 : P — £(/ /i)M o d ? 2 P, 
fc=l i=l 

where M(A® JB(g)C) = AC®B. We also define H = XT 1 -E 2 and S the op- 
erator on C _(X! ,...,X m ) given by HP = E P + ~ X P if P G C (Xi , . . . , X m ) . 
We extend 3 to C(Xl, . . . , X m ) by setting HI = 0. We set, for i = 0, 1, 2 or 
nothing, S; the operator on C(Xi, . . . , X rn ) m such that Si(Pi, . . . , P m ) = 

("iPl j • • • j ^jPm) • 

Lemma 5.1. Por aZZ Z G {1, . . . ,m}, for all P G Co(Xl, . . . , X m ), the fol- 
lowing equalities hold: 

D l T,- 1 P = T 1 - 1 D l P + D l P, 

m 

AHiP = Si ASP + diDMD&P, 

i=l 

Di~ 2 P = E 2 ASP. 

Besides, let Hess(F) : C{X U . . .,X m ) m C(X l5 . . . ,X m ) m 6e groen fry 

m 

Be SS (V)(v) l = Y,d l D l V^v i . 

i=l 

Then, for any (Pi, . . . , P m ) G C(Xi, . . . , X m ) m , with I the identity on C(Xi, 
. . . , X m ) m , the following relation of commutation relation holds: 

DE = (I + Hess(y) + 3)DS. 
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Proof. By linearity, it is sufficient to prove these equalities for a mono- 
mial P = Xi 1 ---Xi p . Moreover, the projection II onto Co(Xi, . . . , X m ) is 
irrelevant in the definition of the operators Sj's since they are followed by 
derivatives: 

D^P = P D t P = (p - l)D t P + D t P = S- 1 DtP + D t P. 
To prove the second equality, write 

D l E l P = D l Y qiDiVqs, 

i,T,P=q 1 X i q 2 

then Di can differentiate qi, q2 or DiV so that 

D{E\P = Y r 2 D i Vr 3 r 1 + Y r 3 r 1 D i Vr 2 

i,Y,P=riXir2Xir3 i,T,P=riX ^2X^3 

+ Y <?4<?2gi93- 

i,T,P=q 1 X i q2,D i V=q 3 X l q i 

The sum of the first two terms gives exactly Si-DjEP and the last one is 

Y qiDiPqs = diDtVtDiZP. 

i,D i V=q 3 X l q i 

Note that if P is a monomial, 

E 2 P = 2 Y M<?i<?3]<?2 + /4<?2]<?i<?3} 

i,'EP=qiXiq2Xiq3 

so that we obtain 

D{~ 2 p = 2 Y vfaWiqsqi 

i,'EP=qiXiq' 1 X i q2X i q3 

+ 2 Y v[qzqi]q'2<& 

i,SP=qiX i q2Xiq' 2 X i q3 

+ 2 Y Mtekkiftj- 

i,'EP=q 1 X i q2X i q3Xiq' 3 

Similar algebra shows that 

S 2 ASP = 2 Y {fi(q 2 )q 3 qi} = D^P 

i,DiXP=q 1 X i q 2 X i q 3 

Finally, the last point we only have to sum the previous equalities for P £ 
Co(Xi, . . . ,X m ) and all I € {1, . . . , m}, 

m 

Dt-Z-ip = (A + 1 A - H 2 A + Hi + E diDtV^DiP 

i=l 

= [(/ + Hess(y ) + S)DP], . □ 
Thus, we can deduce an expression for D o EH -1 . 
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Lemma 5.2. The operator E is a symmetric nonnegative operator in 
L 2 (ii). Let •* be the involution on C(X\, . . . ,X m ) <g> C(X±, . . . ,X m ) defined 
by (A ig) B) 1 = B® A, then for any (P,Q) G C(Xi, . . .,X m ) m , 

in 

Lt(PEQ) = Y, x [d fe Q]*). 

fc=i 

S is i/ras nonnegative in L 2 ( / u) m equipped with the scalar product (P,Q) = 
n(Pi.Qi) ■ + Hess V is a nonnegative operator in the sense that for 
every polynomial Pi, ... , P m , 

m m 

^(Hess(y)P),P* > -(1 - c) £ 

i=l i=l 

Thus, (I + Hessy + H) is symmetric definite positive in L 2 (fi) m and is in- 
vertible. If we consider DSH -1 as a continuous operator from Co(-Xi, . . . ,X m )A 
into L 2 (fi) m , the following rule of commutation holds: 

DSH" 1 = (I + Hess V + H) _1 D. 



Proof. Here, it is easier to come back to the origin of the problem. The 
idea is that the operator H is a projection of the Laplace operator 

-, m N 
k=l i,]=l 

on functions of the matrices. Here, cLfc is a notation and stands for 
In fact, if we take P a polynomial function, 



LP = 1 

N 



m N m N 

-k=li,j=l k=li,j=l 



= E d k Pt(-D k V -X^ + ^I® ft N )M o (d k o d k )P, 

k=l k=l 

with Ajj the matrix with null entries except in where it is equal to 1. 
As a consequence, we deduce from the convergence of £i N toward fi that, for 
all polynomials P, Q, 



lim / ^TV(QLP) di% = -m(QSP). 

TV— >oo J iv 
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But now, by integration by parts, we obtain 
1 

N 

N 



: Ti{QLP)d^i v N 



f 1 

= / ^ E QaA LP )l3,adll 
a,/3=l 

r -, m N 

(20) =-JW^ £ d^Qa^Pp,^ 

k=l i j',a,/3=l 
fc=l i,j,a,0=l 

m . 

= U N ® fa N {dkP x (£> fe Q)') d^, 

k=i J 

which converges as N goes to infinity toward 

m 

® KdkP x (d k QY) = fjt(QEP). 

k=l 

This shows that 5 is symmetric and nonnegative [since if Q = P*, the right- 
hand side of (20) is clearly nonpositive for all N]. Similarly, remark that 

(HessyP) z = £>AnP;- 

i 

Once estimated at a finite matrix, it is easily seen that 

WW)= E (d xi d xl TrV^Pi)^^ 

and so the positivity of Hess is deduced at finite iV from the convexity of V 
which, by definition, is the positivity of the Hessian of Tr(V) in any finite 
dimension. As a consequence, the operator / + Hess(V) + S is invertible on 
C(X\, . . . : X m ) m C (L 2 (/i)) m . We then obtain the commutation relation by 
using the third point of the previous lemma. □ 

This gives us an explicit formula for a 2 . 

Lemma 5.3. For all P, Q in Cq(Xi, . . . ,X m ), for all 1 < k,l < m, the 
following identities hold: 

H ® ^{dkDtP x diD k Q] = /i® ufaDiP x [dkDiQ]% 

m 

C(DP, DQ) = KAP[(/ + Hess V + H)DQ] 4 ) , 

i=l 
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m 

a 2 (EP, Q) = Y^ m(ASPAQ) , 

i=i 

m 

a 2 (P, Q) = £ KAP(I + Hess V + 3)" 1 DiQ). 

i=i 

Proof. An elementary computation shows that, for all polynomials P, 

d k D l p = (d l D k p) t . 

To prove the second equality, recall that 

rn 

C(DP, DQ) = £ (A* ® M[9 fe A^ x ^Z) fc Q] + o 8 k V$(D k P, D t Q))) 

k,l=l 
m 

+ Y / KD k PD k Q). 

k=l 

The third term can be directly written Yli=i /i(A^[DQ]i). For the second 
term, we use the first equality and Lemma 5.2: 

m m 

J2 x diD k Q] =^2fj,(DiPEDiQ). 

k,l=l i=l 

Finally, we only need to check if the two terms in the second derivative of 

V coincide, but this is clear by the trace property: 

m m 

n{dl°d k V%{D h P,D l Q))= KDiPdjDiVftDjQ). 

k,l=l i,j=l 

For the last points we only have to use the commutation rule of Lemma 
5.2 and the previous point: 

a 2 {EP,Q) = C(DSP,DSH" 1 g) 

m 

= /i(ASP[(J + Hess V + S)DSH" 1 Q] i ) 

i=i 
m 

= $>(AEPAQ). 

i=i 

The last point is proved with the same technique. □ 

5.2. Combinatorial interpretation. It was shown in [14] that for small t's 
the limit measure \x has a combinatorial interpretation. More precisely, let 

V = J2i UQi with some monomials qi. Note that in order to have a self-adjoint 
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potential, in the decomposition in monomials, the coefficient of a monomial 
must be the complex conjugate of the coefficient of its adjoint. 

We define a set of colors as the set {1, . . . , m} and associate to each mono- 
mial q = Xi x ■ ■ ■ Xi p a star (i.e., a vertex with some half-edges pointing out 
of it) of p half-edges which are in the clockwise order respectively of color 
ii, «2) ■••)*?■ Besides, we distinguish the first half-edge so that we clearly 
obtain a bijection between monomials and stars. We will say that the star 
is of type q if it comes from a monomial q in that way. Note that a star can 
equivalently be represented by an annulus with ordered colored dots and a 
distinguished dot. 

Given a set of such stars embedded in the sphere, we can construct some 
graphs among them simply by gluing pairwise different half-edges of the 
same color and such that the resulting edges do not cross each other. We 
call a graph obtained in this way a planar graph. Two planar graphs are said 
to be equivalent if there is a homeomorphism of the sphere which fix each 
star and take the first graph on the second. A map is a class of equivalence of 
connected planar graphs for the relation of homomorphism. We now define 



and 



\a r p\ _ u / ma P s with ki stars of type qi 1 

fcl -' fe ^ j-B \ " and one of type P J 

ka (P Q) — if ma P s with h stars of type qi 

ki,...,k n { ,<J) — B | Qne of type p and Qne of type Q 



These quantities are only defined for P and Q monomials, but we im- 
mediately extend them by linearity to arbitrary polynomials P and Q. 
By convention, the star associated to the monomial 1 is empty so that 
M kl> ..„ kn (P,l) = 0. 

In [14], Section 3.2 there is the following relation between the limit mea- 
sure and the enumeration of planar graphs. 

Theorem 5.4. There exists ij > such that, for t S B v , for all polyno- 
mial P, 

MC)= E Il^M^iP). 

k\ , . . .fen % — 1 

We now prove that there is a similar link between the variance o~ 2 (P) 
which appears in our central limit theorem and the generating function of 
theM kl ,..., kn {P,Q)- We define 



n (-t ) ki 
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We shall prove that a 2 (P, Q) and A4(P, Q) satisfy the same kind of induction 
relation. 

Proposition 5.5. For all monomials P,Q and all k, 
M kl ,..., kn (X k P,Q) 

= E E UC%M pl ,... J *{R,Q)M kl -n,..., ka - Pn (S) 

0<Pi<ki P=RX k S i 

+ E E U^ M Plr..,Pn(S,Q)M kl - pi ,..., kn - Pa (R) 

0<p % <kiP=RX k S i 

+ k j M klt ^ kj - 1 ,^ kn {D k VP,Q)+M kl ,^ kn (D k QP) 

0<j<n 

and 

(21) M(X k P, Q) = M((I ® n + n ® I)d k P) - M{D k VP, Q) + n(D k QP). 

Besides, there exists r\ > so that, there exists R < +oo such that for all 
monomials P and Q, all t £ B(0,r]), 

\M(P,Q)\<R dcsP+dcgQ . 

Proof. The proof is very close to that given of Theorem 2.2 in [14] 
which explains the decomposition of planar maps with one root. We look at 
the first half-edge with color k corresponding to in X^P: 

1. The first possibility is that the half-edge is glued to another half-edge of 
P = RX k S. It cuts P in two monomials R and S and it occurs for all 
decomposition of P into P = RX^S which is exactly what does D. Then 
either the component R is linked to Q and to pi stars of type qi for each 
i, this leads to 

U^kiMpu--., Pn (R,Q)M kl - P1 ,... !kn - Pn (S) 

i 

possibilities, or we are in the symmetric case with S linked to Q in place 
of R, 

2. The second case occurs when the half-edge is glued to a star of type qj 
for a given j, then first we have to choose between the kj vertices of this 
type, then we contract the edges arising from this gluing to form a star 
of type DiqjPi; there are 

kjM ku ... :kj - li ... tkn (D k q j P,Q) 

choices. 
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3. The last case is that the half-edge can be glued with the star associated to 
Q = RXiS. We contract this half-edge and obtain a star of type D k QP. 
This leads to 

M kl ,..., kn (D k QP) 

possibilities. 

We can now sum on the k's to obtain the relation on A4. 
Finally, to show the last point of the proposition, we only have to prove 
that there exists A > 0, B > such that, for all /c's, for all monomials P and 

Q, 

M kl ,...,k n (P,Q) < jS^.k, ndegP+degQ 

This follows easily by induction over the degree of P with the previous 
relation on the A4 since we have proved such a control for M. kl kn (Q) i n 
[14]. □ 

We can now relate the variance and the generating function for the enu- 
meration of planar maps with two prescribed vertices. 

Theorem 5.6. Assume (H') with rj small enough. Then, for all polyno- 
mials P, Q, 

a 2 (P,Q)=M(P,Q). 

Proof. First we transform the relation on M. We use (21) with P = 
D k T,R to deduce 

M(ZR,Q) =J2»(DkQD k ZR). 

k 

Let us define A = a 2 — M.. Then according to (5.2) and the previous prop- 
erty, A is compactly supported and for all polynomials P and Q, 

A(EP,Q) = 0. 

Moreover, with M(l, Q) = = cr 2 (l, Q), 

A(1,Q) = 0. 

To conclude, we have to invert one more time the operator H. For a polyno- 
mial P, we take, as in the proof of the central limit theorem, a sequence of 
polynomial S n which goes to S = H~ 1 P in Co(-X"i, ■ • ■ ,X m )A- Then, write 

A(P, Q) = A(E(S n + S-S n ),Q) = A(S(5 - S n ),Q). 

But by continuity of S, S(5 — S n ) goes to for the norm || • ||^. We can 
always assume A > R if 77 is small enough. Moreover, because A is compactly 
supported, A is continuous for || • \\a, and so A(H(5'— S n ),Q) goes to zero 
when n goes to +00. This proves the theorem. □ 
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6. Second order correction to the free energy. We now deduce from the 
central limit theorem the precise asymptotics of N5 (P) and then compute 
the second order correction to the free energy. 

Let 4>o and (j> be the linear forms on Co(-Xi, • • • ,X m ) which are given, if 
P is a monomial by 

m 

(22) mp) =£ £ ^(ftft.ft) 

i=i P=Pl x 1 p 2 x i p 3 

and <j> = (fio o S. 

Proposition 6.1. Assume (H'). T/ien, /or any P in <Co(Xl, . . . 
lim N6 N {P) = 0(H _1 n(P)). 

TV^oo 

Proof. Again, we base our proof on the finite dimensional Schwinger- 
Dyson equation (13) which, after centering, and since we can always assume 
that P £ Co(X\, . . . ,X m ), reads for i G {1, . . . 

N 2 ^((fi N - + AVOP - (/ ® m + J)$J°D = ® 5 N (diP)). 

Taking P = P>jXP and summing over i G {1, . . . , m}, we thus have 



(23) N 2 ^((ft N - M )(HP)) = M # ® ^ (£ ^ o DjXP 
By Corollary 4.8 and Lemma 5.1, we see that 

Jm^f^g^ffflioASP]] =0(P), 



u=l 
?TV. 



which gives the asymptotics of Ao" (HP) for all P. 

To generalize the result to arbitrary P, we proceed as in the proof of 
the full central limit theorem. We take a sequence of polynomials Q n which 
goes to Q = E _1 P when n goes to oo for the norm || • \\a- We denote R n = 
P — EQ n = S(Q — Q n ). Note that as P and Q n are polynomials, then R n is 
also a polynomial. Then we write 

NS N (P) = N5 N (EQ n ) + N5 N (Rn). 

According to Property 3.1, for any monomial P of degree less than eA 2 / 3 , 

\N5 N (P)\<C des{p) . 

So if we take the limit in A, for any monomial P, 

lim sup \N5 N {P)\ <C dc s( p ) 

TV 
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and if P is a polynomial, 

limsup|A^(P)| < \\P\\ C < ||P||a- 
N 

The last inequality comes from the hypothesis (H') which require C < A. 
We now fix n and let N go to infinity, 

limsup|iV5 7V (P -EQ n )\ < \imsup \N S N (Rn)\ < \\R n \\ A - 

N N 

If we now let n go to infinity, the right-hand side term vanishes and we are 
left with 

\imNS N (P)=limlimN8 N (Q n ) = lim0(Q n ). 

It is now sufficient to show that <j) is continuous for the norm || • \\a- 
But it can be checked easily that P — ► Y^=\ d% ° D%P is continuous from 
Co(-Xi, . . . ,X m )A to Cq(Xi, . . . , X m )A-i and a 2 is continuous for || • \\a-i 
due to the technical hypothesis in (H'). This proves that cf) is continuous 
and then can be extended on Co(Xl, . . . ,X m )A- Thus, 

UmN5 N (P) = limc£(Qn) = HQ)- □ 

TV 7i 

This result allows us to estimate the first order correction to the free 
energy. 

Theorem 6.2. Assume (H'), then the following asymptotics hold: 
logZ$ = N 2 F°(V t ) + F^Vt) + o(l), 

with 



F°(V t ) 



J Q f^at^Ylkqij da 



and 



F 1 (V t ) = - £ Kt [rT a \ E kq^J ds, 



with H Q t (resp. (j) a t) the operator H (resp. the linear form 4>) corresponding 
to the potential V a t = aVt with parameters at. 

PROOF. Remark that, for i £ {1, . . . , n}, 

d a log Z\ = -N 2 ^ Vt (± ) 
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so that we can write 



(24) \ogZ^ = N 2 F\V t )- f 

J 



da. 



Since for all a G [0, 1], V a t = aVt is cA 1-convex if Vt is c-convex, Proposition 
6.1 and (24) finish the proof of the theorem since, by Proposition 3.1, all 
the N5^ t (qi) can be bounded independently of N, a E [0, 1] and t S B V:C so 
that the dominated convergence theorem applies. □ 

As for the combinatorial interpretation of the variance, we relate F 1 ^) 
to a generating function of maps. This time, we will consider maps on a 
torus instead of a sphere. Such maps are said to be of genus 1. We define 

■ jl ,ps _ a f maps of genus 1 with k{ stars of type qi \ 
ki,...,k n \ ) — fl I and one of type P J 

and 

fc n = (J{maps with fcj stars of type qi}. 
We also define the generating function 

If P is a monomial, we will denote M.(diP) for X)p=itx i s-M(-^>^>) and we 
extend this notation to all polynomials by linearity. 

Proposition 6.3. For all monomials P and all k, 
Ml^jX.P) 

= E E n^L.., P ^)^pi,..^- P .( s ) 

0<Pi<k t P=RX k S i 

+ e e n^pi,..*.^*!-*!,..^-*.^ 

0<pi<kiP=RX k S i 

+ E + E ^*l,..,fcn(^^ 

0<j<n P=RX k S 

and 

(25) M^XfcP) = 7W X ( (I®ii + fi(g) I)d k P) - M l {D k VP) +M(d k P). 

Besides, for r/ small enough, there exists R < +oo such that, for all mono- 
mials P, all t G -6(0, rf), 

IM^P)! <P de s p . 
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Proof. We proceed as we did for the combinatorial interpretation of 
the variance. We look at the first half-edge corresponding to X k , then two 
cases may occur. 

1. The first possibility is that the half-edge is glued to another half-edge of 
P = RXfrS. It forms a loop starting from P. There are two cases: 

(a) The loop can be retractable. It cuts P in two monomials R and S and 
it occurs for all decomposition of P into P = RX k S which is exactly 
what does D. Then either the component R or the component S is 
of genus 1 and the other component is planar. It produces either 

i 

possibilities or the symmetric formula (where we exchange R and S). 

(b) The loop can also be nontrivial in the fundamental group of the 
surface. Then the surface is cut in two. We are left with a planar 
surface with two fixed stars R and S. This gives 

M ku ..., kn (R,S) 

possibilities. 

2. The second possibility occurs when the half-edge is glued to a half-edge 
of a star of type qj for a given j, then first we have to choose between the 
kj stars of this type, then we contract the edges arising from this gluing 
to form a star of type DiqjPi; this creates 

k j Ml u _.. tk .^.. tkn (D k q j P,Q) 

possibilities. 

We can now sum on the fc's to obtain the relation on Ai 1 . 
Finally, to show that M 1 is compactly supported, we only have to prove 
that there exists A > 0,B > such that, for all fc's, for all monomials P, 

M k,-,k n ( P ) < A J2.k iB de S P 

Another time this follows easily by induction with the previous relation on 
the M 1 (P)'s. □ 

We now give the combinatorial interpretation for the first order correction 
to the free energy. 

Proposition 6.4. Assume (H'). There exists r/ > small enough so 
that, for t 6 B VtC , for all nonconstant monomial P, 

cp(~T l P)=M l {P) 
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and 



n <-t ) ki 

**= e n^*c.,*»- 

fei,...,fe„eN™-{0} i=l *' 

Proof. We use the previous property with P = D k T,P and we sum on 

k: 

M\EP) =M[J2 d k D k EP) =J2<r 2 (dkD k XP) = <f>(P), 

\ k J k 

where we have used the combinatorial interpretation of the variance (The- 
orem 5.6). As M 1 and 4> are continuous for || • \\a when rj is small enough, 
we can apply this to zT 1 P and conclude. 

Finally, for rj sufficiently small, the series is absolutely convergent so that 
we can invert the integral and the sum to obtain 

FHVt) = - f Q M l atl _ atn da 

=[ e En L t (-*i)^ 1 ,.^fe)^ 

J0 fci,..,fen 3 i % - 

= E —i— ... —i— \c 1 1 E II ■ rj (- t j) M l 1 ,...,k J +i,...,k n 



HO 



k i 



e n p M k,...,k j ,...,k„- 

u 

7. Diverging integrals. Physicists often use matrix models in more gen- 
eral settings. We would like to study the case of a potential V for which 
the integral Zy is not convergent. For example, one may wonder if we can 
obtain the generating function for planar triangulation. The issue is that for 
V = iX 3 , Zy is infinite. The idea to give a meaning to this integral is to 
add a cut-off; we define, for L > 0, 

Z$ L = I e —NTr(V(Ai,...,A m )) ^ 

max 

(A)<L 

This allows us to define the probability measure 
Hy AdAi, . . .,dA m ) 



t 



X^(A)<L e - NT r{V(A lr ..,A m )) ^ ^ 
Z V,L 
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In [14], we show that, for all L > Lq for a well chosen Lq, there exists r/ > 
such that for |t| < r], fi N goes almost surely toward the unique solution 
to Schwinger-Dyson's equation (2). This shows that the cut-off does not 
perturb too much the model since the limit does not depend on the choice 
of the cut-off L and keeps the same interpretation than in case of convex 
potentials. The aim of this section is to show that we can also extend the 
central limit theorem to this setting. The key idea is to see this potential as 
a convex potential. We bound the Hessian of 

(26) (f$ t : (A k (ij)) € (R N2 ) m n {A max (A) < L} -> Tr(V(Ax,. . . , A m )) 
uniformly in N: 

n 

Ress ip$ t (A, A) = Y,U Yl Tr(RASAT). 

i=l q i= RXSXT 

Now, using Holder's inequality, 

\Tr(RASAT)\ = \Tr(TRASA)\ < ^Tr ( (TR)A* A (TR)* ) ^Tr (SA* AS* ) 
< ||Ti?||||S||Tr(A,4*), 
which implies that, for {A max (A) < L}, 

II Hess tp$ t || < C|t| 

and C depends only on L. Therefore, we can find e > such that if t E 
B(0, e) n {t| V t = V t *}, for all N, ^ + \ ELi Tr(Xf) is convex on {A max (A) < 

Thus, Vt(A) = Vt(A) + oolA max (A)>L is a convex potential and 
tx mMA) < L e- N ^ V ^=e- N ^A)) 

is log-concave so that most of the step we proved so far can be generalized to 
this case. Indeed, the Brascamp-Lieb and concentration inequalities do not 
require smoothness for the potential V . In fact, we could have included this 
case in all of the previous proofs but they would have been less readable. 
We will only sketch the proof in this generalized case and highlight the main 
differences with the convex case. 

First, we must control the rate of convergence of the measure to its limit. 
The important fact is that up to the choice of t we can obtain bounds 
independent of L. 

Proposition 7.1. There exist nonnegative constants Lo,Mo,C,a such 
that, for L> Lq, we can find rj > such that, for |t| < 77: 
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1. 



2. For all M > M 



n{Xf n ) < limsupJ[ N (Xf n ) < C 2n . 

N 



,,Nf\N /iK A/T\ ^ „-aMN 
^v( A max( A ) > M ) <e 

3. There exists a finite constant sp 1 M such that, for any e > 0, 
M£({|I>) - m N PM \ > e + e% M } n A&) < Jte"^ 



and if P is a monomial of degree d, £pm < NCdM d e 



aMN 



Proof. Since e~ ^^K^C-A)) is log-concave, we can still use the Brascamp- 
Lieb inequalities. The only point to check is that we can still find a lower 
bound for Zy L , but this was already done in [14] using Jensen's inequality: 



JA rT1! ,,x(A)<L 



'A max (A)<L 

II dfi N (Ai 



> ^(A max (A) < L)expf-iV f TV (14(A)) - 



<(A)<L " ' ' V 7V (Amax(A) < L) 

The biggest eigenvalue goes almost surely to 2 and 

1 

{ (A)<L 

is bounded by /i 7V (V t V t *) 1//2 which goes to cr^VtVj*) 1 / 2 < +°o according to 
[23]. Thus, if L > 2, > e~ dN2 for a finite constant d. Thus, we can prove 
the property as in Section 2. The proof of the two last points do not differ 
from the convex case. □ 

The idea, once we have an a priori control on the radius C of the support 
independently of L, is that we can use it to approximate any polynomial by a 
compactly supported function with support in [— L, L]. We choose L> Lq = 
max(Mo,C) and define for Lq < R < L, <f>R the piecewise affine function 
such that, for |x| < R, 4>r(x) = x and (f>R has a compact support strictly 
inside [— L,L]. Then we can approximate any polynomial -P(X) by hu = 
P(<pn(Xi), . . . , (j>n(X m )). The main improvement in the replacement of P 
by ha is that hn satisfies the finite Schwinger-Dyson's equation (13). 

Proposition 7.2. If L is bigger than some L > 0, and e > 0, there 
exist C,t],Mq such that, for M > Mq, |t| < rj for all polynomial P of degree 
d<eN 2 / 3 , 

\6 N (P)\<C^. 
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Proof. In order to prove the analogue in the convex case (Property 3.1), 
we use the finite Schwinger-Dyson's equation which is not always satisfied 
in this case. In fact, it is only satisfied for compactly supported function 
h with support in [—L,L], since for such h we can make the infinitesimal 
change of variable. For a polynomial P, 

${(L N [{X t + DiV)P\) - fi$(fi N ®fi N (diP)) 

= $ (fr N KXi + DiV)(P - h R )}) - ^(fi N fi N (di(P - h R ))). 

Therefore, since \x satisfies the Schwinger-Dyson equation, we get that, for 
all polynomial P, 

(27) 6 N (XiP) = -S N {DiVP) + 5 N ®]l N (d l P)+ f i(g)5 N (d l P)+r(N 1 P), 
with 

r(N,P):=N- 1 ^(f ^f(diP)) 

+ N(^(fl N [(X i + D i V)(P-hR)}) 
-^(^®^( dl (P-h R )))). 

Thus, the only difference with the convex case is the term N{${ji N [{X i + 
DiV)(P - h R )]) - Hv(» N ® fi N (9i(P ~ hit)))) but, since on A^, P(A) = 
h R (A) and R > M, this term decreases exponentially fast and this allows to 
finish the proof exactly as in the proof of Proposition 3.1. □ 

Since the main tools are available, we next turn to the proof of the cen- 
tral limit theorem. Here we have to be careful since the technique of the 
"infinitesimal change of variable" is no longer true in its full generality. But 
it still holds if we restrict ourseves to compactly supported functional, thus, 
we immediately obtain a weaker version of Lemma 4.1: 

Lemma 7.3. If L is bigger than some Lq > 0, there exists r/ such that, 
for |t| < r\ if hi, ... , h m are compactly supported with support in] — L,L[ and 
self-adjoint, the random variable 

m 

Y N (h u . . . , hm) = NY,{^ N ® P^idkhk) ~ V N [(X k + D k V)h k ]} 

k=l 

converges in law toward a real centered Gaussian variable with variance 



m m 

C(hi,...,h m )= {^®n[dkhi x dih k ] + fi(di o d k V$(h k ,hi))) + Y fi(h 2 k ). 

k,l=t k=l 
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The last step is to show that even if we do not have the result for all Stieljes 
functions, it is sufficient to approach polynomials by compactly supported 
function with support inside ] — L,L[. We will again use the fact that the 
limit measure has a support bounded independently of L. Using this idea, 
we prove a result similar to Lemma 4.2. 

Lemma 7.4. If L is bigger than some Lq > 0, there exists rj such that 
for |t| < rj, if P\,. . . , P m are in C(Xi, . . . , X m ) sa , then the variable 

m 

Y N (P l ,...,P m ) = Nj2 \fi N ® t^idkPk) ~ p?[{X k + D k V)P k }\ 

k=l 

converges in law toward a real centered Gaussian variable with variance 

m m 

C(P U ..., P m ) = £>® v[d k Pi x 8iP k ] + nidt o d k V#(P k , P^)) + Y, KPk)- 

k,l=l k=l 

Proof. First choose L > Lq = max(Mo, C) and for Lq < R < L, approx- 
imate the polynomials Pi(X.) by h l R = Pi((j)R(Xi), . . . , (f>R(X m )). Then since 
C bounds the support of /x, observe that C(P\, . . . , P m ) = C{h l R , . . . , h R ) and 
we only have to prove that 

Y N (P 1 ,...,P m )-Y N (h 1 R ,...,h^) 

goes in law to when N goes to infinity. But, we have the inequality 

P(\Y N (P 1: . . .,P m ) - Y N (h R , . . . ,/#)| > e) < P(A max (A) > R) 

and the right-hand side goes exponentially fast to 0. □ 

The other results can be proved as in the convex case with only minor 
modifications. Following the same way than in the convex case, this allows 
us to prove the theorem: 

Theorem 7.5. If L is bigger than some Lq > 0, there exists rj such 
that, for |t| < r), for all P in Co(-Xi, • • • ,X m ), 5 N (P) converges in law to a 
Gaussian variable with variance 

a 2 (P) := Ci-^P) = CiD^E^P, . . . , D^E^P). 

Besides, the convergence in moments occurs and the covariance keeps its 
combinatorial interpretation, allowing us to enumerate a larger variety of 
graphs. 

Finally, applying the same strategy than in the convex case, we are able 
to prove the convergence of the free energy. 
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Theorem 7.6. For L is bigger than some Lq > 0, there exists r/ such 
that, for |t| < rj, the following asymptotics hold: 

log Z$ uL = N 2 F°(V t ) + F\V t ) + o(l), 

with 

fel,...,fc„GN-{0}i=l l ' 

and 

fcl,...,fc„GN-{0}i=l l ' 
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